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FORMULAE AND DIAGRAMS FOR THE CORRECTION OF THE SYSTEMATIC 
ERRORS OF THE HELSINGFORS ZONE OF THE 
ASTROGRAPHIC CATALOGUE 


By HEINRICH EICHHORN 
Georgetown College Observatory, Washington, D. C. 
Received September 11, 1958 


Abstract. The star positions in the Helsingfors zone of the Astrographic Catalogue are affected by systematic errors 
(Ag, An) that depend on the magnitude and the position of the stellar image on the plate. These errors will, in extreme 
cases, produce systematic differences up to 150 between the positions of a star of sixth and one of tenth magnitudes. 
This paper gives formulae for the corresponding systematic corrections in units of o‘o1 to the standard coordinates 
‘t and 7, in the form —AE = (m — 10™0)(XM — 4.4) + X1, and —Ayn = (m — 1070) YM + Y1, where m is the directly 
estimated magnitude denoted by ‘“‘Gr’’ in the Catalogue. Diagrams are provided for the rapid evaluation of XM, X1, 
YM and Y1 as functions of the star’s measured positions x and y in millimeters. The application of these corrections, 
which will increase the weight of the positions in extreme cases by a factor of 50, is especially important when star posi- 


tions are computed from one plate only. 


I. INTRODUCTION, GENERAL REMARKS ON 
THE COMA EFFECT 


The Helsingfors Zone of the Astrographic 
Catalogue (ACHels) is one of the most accurate 
zones. In their investigations on accuracy and 
systematic errors of the positions, Furuhjelm 
(1906) and Groenstrand (1937) found the scale 
of the plates to depend on the magnitude. A 
systematic error of this type, probably due to a 
slight over- or undercorrection of the objective 
lens for coma, was also found by Schlesinger and 
Hudson (1914) on plates taken with a wide angle 
doublet and by Schlesinger (1925) on the plates 
taken for the Yale zones. Eichhorn (1956) found 
the same effect to be rather large, and possibly 
variable with seeing, on plates taken with the 
26-inch refractor of the Leander: McGormick 
Observatory. A strong coma effect also exists in 
the Hyderabad zone of the AC (Eichhorn, 1957) 
where, in extreme cases, the systematic errors 
amount to several seconds of arc. 


2. SYSTEMATIC ERRORS IN PHOTOGRAPHIC 
STAR POSITIONS 


The deviations from a strictly gnomonic rela- 
tionship between a star’s spherical coordinates 
and its standard coordinates as obtained from 
the measured coordinates of its image on a photo- 
graphic plate will be functions of its coordinates 
on the plate, its magnitude and color index, the 
optical properties of the objective, the sensitivity 
function of the film, deviations of the film from 
a plane and from the focal surface, guiding and 
the atmospheric conditions during the exposure. 


A complete theoretical evaluation of all these 
effects would be very involved if not forbidding, 
and only a more or less empirical approach seems 
possible at present. 

Disregarding the influence of the atmospheric 
conditions, we can in principle write the correc- 
tions for the reduction to gnomonic projection 
in the following schematic form: 


AE = U(m,c) + Vim, c, x,y) + W(x, y) 


and an analogous expression for An; here x and y 
are the rectangular coordinates of a star on the 
plate, m is the magnitude and c a measure of its 
color. To determine expressions for Aé and: Ay 
we set up trial formulae for U, V and W with 
enough parameters to keep them sufficiently 
flexible, and then determine these parameters 
empirically. In the absence of more specific. in- 
formation we assume that the systematic errors 
(i.e, U, V and W) are numerically the same for 
all the plates taken with the same instrument. 

There are two ways of establishing observation — 
equations for the parameters in the functions U, 
V and W: 


a) Comparison of stellar positions from the 
plates with positions from a standard catalogue 
(Catalogue comparisons). These are capable of 
providing the parameters of functions of any of 
the types U, V or W. 

b) Comparisons of stellar positions as obtained 
from partially overlapping plates. Let x1, y, be 
the coordinates of a star on the one plate, x2, ye 
the coordinates of the same star on the other 
plate. As m and c are characteristic for the star 
regardless of its position on the plate, U cannot 
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be determined from material of this kind. Fur- 
ther restrictions may occur as consequences of 
the relative positions of the stars on two plates. 
If, for example, |x1| & |xe|, |1| & |yel, cer- 
tain parts of Vand W will remain undetermined, 
for instance those which depend on x? or y?. 


In the case of intercomparison of plates the 
observation equations will have the form 


6& a Agi ae Ag, = Vim, C, X1, V1) 
3 V(m, €, X2, V2) + W (x1, 1) Fes W (x2, 2), 


where the 66, i.e., the differences between the 
corrections, are the directly observed quantities. 

Before evaluating the expressions for the cor- 
rections we shall give an outline of the approaches 
made. Matrix notation is employed wherever 
possible in order to save space. Thus, if the 
corrections are given by 


—AgE = ax + a,x? + ary + --- 
—An = bx + bix? + boxy + --- 


- we write this in matrix form 


Ea) as 
where 
@ ay a2 
4=(5 bbs ) 
and 
x 
x2 


In order to save space we shall indicate column 
vectors like S by their transposes S’, writing 


S’ = (x x? xy ---), 


‘In the case of the Helsingfors zone the avail- 
able data on colors were too scarce to evaluate 
the dependence of the errors on color. We as- 
sumed the errors to be polynomials in the coor- 
dinates, the distance from the plate center 
r = (x? + y*)# and the magnitude. Magnitudes 
m were used after subtracting from the catalogue 
magnitude m, a zero-point magnitude mo, thus 
m= Mm, — mM. m, is always the directly esti- 
mated magnitude found in the catalogues under 
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“Gr.,”’ never the magnitude found under ‘Gr. 
déf.”’ 

Several types of formulae were tried. First we 
used a polynomial in x?, xy, mx, m, I, x with 
mo = 10.0 mag., that is, a coefficient matrix P 
with the transpose 


P! = (xy mxm ix). 


Secondly we tried another coefficient matrix, Q, 
with the transpose 


Q’ = (rmx rmy rm xym tr xy x y 1) 


with my = 8.5 mag. For the final solution, from 
which the correction diagrams were drawn, we 
adopted the coefficient matrix R with the trans- 
pose 


R’= 
(I mx? yxy x?m ym xym xy rXxM rym xm ym), 


and used again my = 8.5 mag. 
The results of the solution can be written less 
clumsily if we use a matrix R,, defined by 


Ry = (1x2 9 xy rx ry x y) 


and then write the solution in the form 


ae! = (mA + B)-R,, 
where A and B are the matrices of the coefficients 
which have to be found empirically from the 
comparison material described below in 3. As no 
terms rx and ry occur in R, the corresponding 
terms in B have to be assumed identically zero. 
By performing the matrix multiplications the 
formulae for the corrections can be brought into 
the following form: 


—AE = m[XM(x, y)] + X1(x, 9) 


From the final solution described below diagrams 
for the functions XM, YM, X1 and Y1 have been 
drawn which depend only on x and y, after 
adapting the formulae to the zero magnitude 
my = 10.0 mag. From the diagrams we can read 
the values of these functions for any pair of x 
and y. This makes it very easy to compute the 
corrections. 


Il 


3. THE MATERIAL 


a) Furuhjelm differences. Differences 5 (unit 
0001) and dn (unit 0’01) of the systematic errors 
Aé and Ay in nine selected points of overlapping 
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quarter plates are given by Furuhjelm (1906) 
for the following groups: field stars (étoiles ordi- 
naires; mean magnitude 10.2 mag.), comparison 
stars (étoiles de rattachement; mean magnitude 
9.7 mag.) and reference stars (étoiles de com- 
paraison ; mean magnitude 8.4 mag.), taken from 
Vol. IV of the ACHels. These error differences, 
obtained by taking means of a total of 4045 com- 
parisons each in a@ and 6, were tabulated against 
the coordinates of one of the two quarter plates 
used for comparison. The coordinates of the cor- 
responding point on the other quarter plate can 
be reconstructed from additional data contained 
in Furuhjelm’s paper. We used the differences 
for the field stars from Furuhjelm’s table on 
p. 8; the analogous data for reference and com- 
parison stars can be reconstructed from the 
equations on pp. 19 and 22. 

b) Groenstrand differences. Groenstrand tabu- 
lates for three magnitude groups (7.2, 8.5 and 
9.6 mags.) means of the differences 5£ (unit of00!) 
and én (unit o’01) of the systematic errors in 
the regions common to plates that overlap on 
their edges and in the centers of overlapping 
plate quarters, against the coordinates on one of 
the plates of the pair (Groenstrand table II, p. 
20). The coordinates of the star on the other 
plate can be calculated. These figures are based 
all in all on 2532 individual comparisons in a and 
the same number in 6 from the ACHels Volumes 
I, III, IV and VIII. 

c) Groenstrand comparisons. Means of system- 
atic errors Aé (unit o*001) and Ay (unit o%or) 
depending on the position on the plate are given 
for four magnitude groups (6.3, 7.5, 8.5 and 9.6 
mags.) and nine normal positions on the plate 
each in a and 6 (Groenstrand table VIII, p. 36) 
based on 4056 comparisons of ACHels positions 
with the meridian catalogue by Morin and Kon- 
dratjeff (1933), which was observed at the re- 

. quest of the Helsingfors observatory for the pur- 
pose of determining the systematic errors of the 
ACHels. The normal positions lie practically on 
the plate diagonals, four of them about halfway 
from the center, four near the corners and the 
ninth at the center of the plate. The systematic 
errors in all nine points were reduced to the:plate 
origin by subtracting the systematic errers at 
the origin from them, so that the errors at the 
origin are by definition equal to zero. No:tse of 
this latter information was made in our compu- 
tations as this circumstance was unfortunately 
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overlooked, mainly because the plate center com- 
parison is not included in Groenstrand’s table 
VIII. A diagram, drawn essentially from this 
material (p. 38) shows clearly that coma causes 
the most important of the systematic errors. 

d) Furuhjelm edge differences. Differences of 
the systematic errors on the edges of overlapping 
plates tabulated against the coordinates on one 
of the plates are also given by Furuhjelm (pp. 
24 and 25). Because of the lack of magnitude 
data for this group, which includes stars of all 
magnitudes, it could not be used in this inves- 
tigation. 

Furuhjelm and after him Groenstrand investi- 
gated their material for errors of the type 


ax+by+cx?+dxy+ey+fxr-+ Aa cos 6=0 
axtbytenx?+diwytey+fyrt+Ab=0. 


Furuhjelm’s analysis, in which cubic distortion 
was not considered, indicated the presence of tilt 
of the plates, absence of quadratic distortion 
(f, f1), and a coma effect. It also suggests that 
the theory of ideal gnomonic projection is in- 
sufficient for the description of all systematic 
effects. The computations of Groenstrand con- 
firm essentially Furuhjelm’s results on the coma 
effect and the optical distortion and make him 
well aware of the fact that the inner systematic 
errors of the positions depend strongly on the 
magnitude, but he declares it hopeless to try to 
remove them. 

The rest of this paper will show that Groen- 
strand’s opinion cannot be upheld. It may be 
that the systematic errors vary to a certain de- 
gree from plate to plate, so that any general 
formula will still leave small systematic errors; 
but the greater part of the systematic errors will 
certainly be removed by a general formula for 
the whole zone and the positions will be ae 
stantially improved thereby. 


4. SOLUTIONS 


All solutions were made by the method of least 
squares. The scheme suggested by Vasilevskis 
(1940) for the solution of normal equations was 
found to be extremely efficient and useful. The 
numerical values for the different coefficient mat- 
rices found, together with the mean errors of 
their elements, are given at the end of this section. 

First a pilot solution for the corrections —Aé 
only was tried with the material from the Groen- 
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strand meridian comparisons in the form — Ag = 
C-P where P is the matrix defined above, and 
C the corresponding matrix of the coefficients, 
the numerical value of which is given below. 
Aé is in units of 0001. The weights were chosen 
to three figures according to the number of indi- 
vidual comparisons that were incorporated in 
any particular normal point, and also considering 
the diagrams given by Groenstrand on pp. 29 
and 30 of his paper. The [pvv ] was reduced from 
3314 to 1592; the mean error of one unit weight 
comparison is 0°078. This rather high value is 
explained by the fact that the comparisons were 
made with 20-year-old meridian positions, mostly 
without applying proper motions, and represents 
rather the error of the meridian positions than 
‘that of the photographic material. Therefore the 
weights of the normal values in the uncombined 
solutions were determined by the number of indi- 
vidual comparisons incorporated in one particu- 
lar normal point. It can be seen that the weights 
would scarcely be changed if the Groenstrand 
weight diagrams were considered. The values of 
Aé, originally in units of o°001, were converted 
to units of 0’01 by multiplication with 1.097. 
Separate solutions were made as follows: 


a) The Groenstrand meridian comparisons 
(equations for 32 normal points) were analyzed 
eed = D.Q. This 
brought the [pvv] down in x from 5543 to 1263, 
in y from 10,453 to 1847, corresponding to a mean 
error of a unit comparison of 07741 in x and 
07896 in y. This material cannot be analyzed 
using the term matrix R, because |x| & |y| for 
all normal points and therefore x? and y* cannot 
be separated. 

b) The Groenstrand differences (78 normal 
points and equations) were analyzed according 
to ea = E,-R = (mA, + B,)-R; and the 
[pov ] in x was reduced from 1408 to 306 and in 
y from 1099 to 231, corresponding to mean errors 
of 066 in x and 0"57 in y. Because of the rela- 
tive situation of the normal points certain un- 
knowns will not occur in the normal equations, 
they are marked by * in the matrices A; and B, 
below. Terms in rx and ry were not included in 
the original matrix R, therefore the value zero 
must be assigned to their coefficients in B,. They 
are the zeros in B, for which no error is quoted. 


according to the formula ( 
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c) The Furuhjelm differences, given in 100 
normal points, were analyzed according to the 
same formula and gave the coefficient matrices 
Ex or Az and By respectively. Elements * and 
zeros without errors mean the same as above in 
b). The [pvv] is reduced from 796 to 209 in x, 
and from 673 to 152 in y, corresponding to mean 
errors for a unit-weight comparison of 0749 in x 
and 0742 in y. 


The agreement within the individual solutions 
is, within the errors, satisfactory. Introducing 
into the final formulae a coma term, quadratic 
in the magnitude m like am?x, might improve the 
aecuracy slightly as the coma seems less pro- 
nounced for the faint stars than for the bright 
stars, but this would probably strain the material 
too much. Furthermore one has to consider that 
the linear magnitude equations are based on the 
comparisons with the Poulkovo meridian cata- 
logue only, which, according to Gyllenberg (1948), 
requires a magnitude correction in right ascen- 
sion of —4.4(m — 4™o) in units of o%01 to reduce 
it to the G.C. 

a ae 
An 


For the final solution in the form (2 
(mA + B)- R all the normal equations were 
added, after those from the Groenstrand com- 
parisons were divided by three and those from 
the Groenstrand differences divided by two, in 


order to reduce the two groups of Groenstrand 


material to the mean error of unit weight of one 
Furuhjelm difference. Thus the weight of the 
10653 comparisons, in x and y each for all three 
groups, is the equivalent of 7166 comparisons of 
unit weight, again in x and y each. The system 
ef normal equations is given in Table I. The 
figures in parentheses are the powers of 10 by 
which the numbers before are multiplied. 

The [pvv] are reduced in x from 1658 to 541, 
and in y from 1571 to 435, corresponding to mean 
errors of a unit weight comparison of 0752 in x 
and 0%47 in y. These are, as expected, somewhat 
bigger than the mean errors of unit weight before 
we put all the normal equations together. The 
fact that the increase is so slight shows that all 
the material is consistent. One should note that 
the terms in xm and rxm on the one hand, and 
ym and rym on the other hand, are very hard to 
separate, so that one cannot really say that the 
separation of the coma in a linear and a quad- 
ratic term is significantly established. 
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TABLE I. THE COEFFICIENTS OF THE NORMAL EQUATIONS 


+135 (1) 1 —40 (1) m —5s9 (4) x4 —62 (4) 98 +4 (4) xy +192 (4) x*m +205 (4) 9°m a 

+117 % +164 . +175 (4 —3 (4 —58 (4) —62 (4) +4 

+2162 (7) +384 (7 +43 (7 +1079 iM —223 (7 +39 

+3043 (7 +46. (7 243507 +1622 (7 +39 

+2970 (7 +40 (7) +41 (7 +1134 

+1448 (7) +255 (7) +8 

1928 (7) +6 

+2468 

—1I5 (2) 9 —5 (4) rxm -—o 4 rym —I0 (2) zm +o (3 ym = 
+0 (2 —6 (4) +3 (4 2 +4 (2) 
—33 (5 —57 (7) —13 (7 —108 (5 —41 (5) 
—39 (5. —2 (7 =4 (7 +12 (5, —12 (5) 
+arr (5 +130 (7, +179 (7 +154 (5 +377 (5) 
—54 (S —39 (7 =13_(7 —64 (5 —33 (5 
—31 (5 +12 (7 —20 (7 +27 (5 —40 (5 
+118 (5 +55 (7 +05 (7 +64 (5 +207 (5 
+108 (3 +8892 (5 +142 (5, +17537 (3 +241 (3) 
+21838 (3 +151 (5 +9399 (5. +241 (3) +18541 (3) 
+0984 (7 +109 a +18140 (5) +193 (5) 
+0621 (7 +179 5} +17689 (5) 
+34206° (3 . +318 2 
+33685 (3 
(pvv] = 
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Table II contains the solutions for the afore-described coefficient matrices with the mean errors 


of the elements: 


_ | —0.00069 
ce [ +167 
—0.00363 0.00000 —-0.00344 
+55 £147 
D = 
0.00000 —0.00249 +0.00564 
+65 +178 
* * 
* * 
* * 
B, = | 
* * 
~! —0.00268 
+96 
A, = 
bY +0 .00236 
+82 
* 0.00022 
ae 
pre 73 
* +0 .00049 
; +63 
+1.9 —0.00297 
1.9 +82 
A= 
+2.8 ++0.00184 
1.7 +74 
—2.7 —0.00014 
+2.0 +66 
B= : 
+3.2  +0.00093 
1.8 60 


TABLE II. SOLUTIONS FOR THE COEFFICIENT MATRICES 


—0.00014 —0.169 
£93 £45 
+0.00082 +5.7 
77 $2.5 
—0.00108 —0.8 
04 +3.0 
‘s +0.00107 
+41 
* —0.00202 
+36 
= +0.00193 
£39 
as —0.00083 
\ £33 
+0.00396 -+0.00182 
+81 #125 
+0.00259 —0.00227 
69 +107 
—0.00074 +0.00175 
+61 +83 
+0.00168 —0.00000 
+52 +75 
+0.00353 +0.00I1I19 
+76 +36 
+0.00208 —0.00198 
+69 £33 
—0.00052 +0.00182 
£59 £33 
+0.00206 —0.00056 
+54 +30 


+1.4 +2.0 +0.19 
£9 


1.5 +3.5 
—0.00236 +0 
£137 
+0.00469 —O 
+164 
—0.0064 —Oo 

+17 
+0.0003 —o 
+15 
0.0000 fe) 
0.0000 Co) 
+0.00091 —O.' 
+164 
—0.00439 +0 
141 
0.00000 
0.00000 Oo 
—0.00160 —Oo. 
+87 
+0.00045 —O. 
+79 
0.00000 Oz 
0.00000 oO. 


.00043 


72 


.00233 


+87 


-OOIO 


£20 


-0037 


+18 


-0000 


-0000 


-00000 , 


+1.8 
+2.3 


+2.3 
+2.8 
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5. THE DIAGRAMS AND THEIR USE for the magnitude equation of the Poulkovo cata- 
logue as quoted by Gyllenberg (1948) ; m, is the 
preliminary magnitude ‘‘Gr.”’ from the catalogue. 
We changed the above-given formulae from the 
— At=(m,—10™0) -[XM (x, y)—4.4]+X1(x, y) zero point magnitude 8.5 which was better for 
Lin ee the computations, to the new zero-point magni- 
ee) Ms, y) FY 1x, 9). tude mo = 10.0, which permits an easier appli- 
The —4.4 in the formula for —Aé takes account cation of the formula. Thus we have: 


The final formulae for the corrections —Aé 
and — An can be written 


XM = +1.9 — 0.00297-x? + 0.00353:7?, + 0.00119:xy — 0.00160-7x — 0.00210:ry — 0.070°% +.0.106-¥ 

XI = +0.2 — 0.00459:x? + 0.00477:y? + 0.00361 xy — 0.00240:7x — 0.00315:77 + 0.035°% + 0.220°¥ 

| YM = +2.8 + 0.00184-x? + 0.00208-y? — 0.00198-xy + 0.00045 :7x — 0.00277:ry — 0.060:% + 0.0II-y 
Yi = +7.4 + 0.00369-x? + 0:00518-y? — 0.00353:xy + 0.00068 rx — 0.00416:ry — 0.005°% + 0.146°y. 


If electronic computers are available these func- tation of these functions for a large number of 
tions are éasily evaluated. However, the compu-_ pairs (x,y) of a desk computer would be very 


+30 8©6+20 +10 fe) 10 -20 -%0 40 -50 -60 


TN JSS Sime 4 \5 - 


Darr reee 


MMEKAIKEL}. : Se 


Og d CE = = at a 
+40 +30 +20—s +O 0 -10 -20 30 40 30 60 


Figure 1. X M: coefficients for evaluating corrections Af. 
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tedious. Therefore we give diagrams which show 
curves X M = 2, where z is any integer that the 
functions X M etc. assume over the field of the 
plate, ie., —65 S (x,y) S 65. 

To get these curves we have to solve equations 
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of the type 
2 = a + ax + dey + asx” + aay” + asxy 

+ aexVx? + y? + aryVx? + 9 


for x or y. Assume we solve for x, then we have 


x = a + ay + aavx? + y? + Vaux? + asy? + ogy (x? + 92)# + a(x? + 9%)! + sy + a — aoe 


which can be used as an iteration formula for x. 
This formula converged rapidly for practically 
all XM and X1 and for a good part of YM. For 
the rest of YM and for Y1 Newton’s iteration 
formula %ni1 = Xn — f(%n)/f’ (xn), properly 


ae 


+40 +30 +10 


adapted for our purposes, was used. The com- 
putations for the diagrams were performed on a 
Burroughs Eto2 electronic computer at George- 
town College Observatory. 

The diagrams give the functions XM etc. in 


-10 -20 


Figure 2. X1: coefficients for evaluating corrections Aé. 
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units of 0701 in a matter of seconds. The im- 
proved standard coordinates are obtained by 
adding the corrections as found from the dia- 
grams to the standard coordinates that are ob- 
tained from the formulae given at the head of 
each plate. Sets of enlarged diagrams with a 
superimposed coordinate grid may be obtained 
from the Georgetown College Observatory. 
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Sere Sd his) 


Figure 3. YM: coefficients for evaluating corrections An. 
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+50 +40 +30 +20 +10 
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THE CALCULATION OF LITERAL EXPANSIONS 


By PAUL HERGET anp PETER MUSEN 
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Abstract. This paper describes a method of using the IBM 650 magnetic drum electronic calculator for making literal 
developments of functions in terms of given parameters. As examples we give series in powers of e obtained for Bessel 
functions, J;(ie), as well as for coefficients arising in expansions of selected functions in elliptic motion expressed as 


Fourier series in the mean anomaly. 


There is a type of problem arising frequently 
in celestial mechanics which requires that a given 
function be expanded literally in terms of certain 
parameters. A simple example is the develop- 
‘ment of the Bessel functions, J;(ée), as power 


Cos 
N +10" x42 X02 x3° X47 5° Xe! X17! sin 


where N is an eight-decimal place number; each 
letter a to n is an integer between —49 and +49; 
and the x’s and A’s may be assigned any literal 
significance desired, so long as it remains fixed 
during any one application. Thus far, our cal- 
culating program can add, subtract, or multiply 
any two terms of this general form. Subse- 
quently, we expect to include differentiation, 


OO IDEN aa bb cc dd ee ff gg 


The sign of the third word is the sign of N, and 
the sign of the fourth word is minus for sine 
terms. The IDENtification is any assigned series 
number from 0000 to 9999. The first two columns 
give the ‘‘order”’ of the term. This is the weighted 
sum of the exponents a to g. The weighting 
factors may be only small integers or zero. 

No attempt is made to accumulate like terms 
within the electronic computer. When two series 
are to be multiplied together, the shorter one 
(limited to 400 terms) is stored on the drum 
first, and then the second series is read from the 
card hopper, one term at a time. The combina- 
tion of two terms, one from each series, produces 
one or two resulting terms, depending upon 
whether at least one of the two is an algebraic 
term or not. A limit may be preset in the calcu- 
lator, and then terms of an order higher than 
this limit will not be punched in the results. 

Furthermore, the limit may be raised when 
certain values of # to / are zero. This is intended 
to serve as a protection against small divisors 
in subsequent integrations. When two series are 


series in e. Other examples will be found in Cay- 

ley’s ‘‘tables.”” We have devised a method for 

calculating such expressions by means of the 

IBM 650 magnetic drum electronic calculator. 
Consider a single term of the form 


(tA, + ide + jAs + kAa + IAs), 


integration, and such other special operations as 
may be needed. It will be noted that algebraic 
power series expansions are included in the above 
general form, if one has cosine terms with h to / 
set equal to zero. 

Each term is carried as a separate IBM ‘“‘half- 
card’’ with the following format for the 40 
columns: 


.NNNNNNNNan hh ti jj kb Ul. 


to be added, they must first be sorted on all the 
arguments a to /. 

The most obvious defect of this whole pro- 
cedure arises from the limitations of floating 
point arithmetic. There are rounding errors of a 
unit in the eighth significant figure, which are 
inevitable, and these are subsequently magnified 
when multiplied by numbers considerably greater 
than unity. We have pondered the prospect of 
computing with only rational fractions, but this 
poses other problems, especially the computing 
time involved in testing every term for the re- 
duction of common factors in numerator and 
denominator. We have decided for the imme- 
diate future to proceed with the present system. 
The experience will make a decision at a later 
time a wiser one. 

In.anticipation of our own future requirements 
for the development of the disturbing function, 
we began by developing the Bessel functions 
Ji-1(te), Ji(te), Jiz1(te) for « = 0 to 20. We re- 
quired two packs of cards containing 1/k! and a 
pack containing (ze/2)* and — (te/2)?. 


12 


Se al 


+ ¢ 16000000 
- «25000000 
+ ¢ 15625000 
643402780 
©67816845 
067816845 
°47095033 
e24028077 
©93859670 
© 28969033 
© 72422583 


fbb ib t + 1 


JP Gee 


+ 616276042 
= 020345052 
+ 010596381 
= 631536851 
+ 061595410 
= 085549176 
+ 689113730 
= 072332574 
+ 047091520 


Je 289 


+ ©4) 704179 
= 093834402 
+ 095007338 
= 058299954 
+ 024595293 
— © 16623802 
+ 018471810 


pla obra Gc} 


0118674718 
© 38592835 
«58233831 
054675098 
© 36094108 


+i etre 


Lom 7 


+ © 35488138 
— © 15082458 
+ 030269657 


Siar Beno 


»50000000 
© 62500000 
© 26041668 
©54253476 
¢67816840 
656514036 
© 33639306 
¢ 15017548 
052144268 
© 14484517 
© 32919360 


+r er eriritis¢ 


51 
50 
49 
47 
45 
43 
41 
39 
36 
34 
31 


THE ASTRONOMICAL JOURNAL _ 


BESSEL 
es ase a Re ee 


+e10000000 
—~50000000 
4083333335 
= 069444446 
4034722222 
~011574074 
#027557320 
~ 249209498 
+ 068346524 
= 075940583 
+ 69036893 


So MEG 


«20250000 
=e 30375000 
+219526786 
=o 73225447 
+ 018306360 
—+32951451 
$244933798 
= 648143354 
$041662520 

TRC Sah 
+ 053822689 
= 613455722 
$015290593 
~ 010618468 
+ 051050328 
= 018232261 
050645166 


J 13 14 


215559449 
= 054458072 
+2e688948161 
=-090801272 
+ 065430328 


dfs) WG ee =) 


+ o46887183 
=~221099234 
+244974681 


th Ue See 


+050000000 
=~ «¢ 16666667 
+ 220833334 
—«13888889 
+257870371 
=¢ 16534392 
+06 34446649 
=054677221 
+068346524 
=269036893 


FUNCTIONS 


Beit tyes 


+211250000 
284375002 
+$023730469 
= 035595704 
4633370972 
=021452769 
210055985 


= 635914232: 


+e10100878 


= e22956539 | 


al jee 


025531467 
=~ 044680067 
4234208178 
=015520376 
+047531153 
=210586484 
4018011727 
= 324246554 


JALO. Ll 


4069801355 
= 019195373 
+024194169 
= 018765990 
+¢10136986 
= 040885841 


J 14 15 


40204386516 
==9 76644431 
+e13472654 
=e 14859545 


J) Lem 


+062041727 
~ 029469820 


Ji eS eee) 


4056250001 
= 031640625 
+ 671191406 
= 688989260 
+ 071509227 
- 040223935 


016759974 


= 053871348 
+ 013773924 
= © 28695675 


52 
53 


50 
50 
49 
48 
47 
46 
45 
43 
42 
40 


JF 3 Ee 


+o 13333334 
—- 0133333346 
+053333334 
=e 11851852 
+e 16931218 
-0 16931217 
+¢12541642 
-0 71666525 
+¢32575693 
=e 12065072 


De Tene 


+032507937 
—065015872 
+¢57791889 
—« 30822341 
+011208123 
—~ 29888327 
+061309395 
-—« 10009697 


Jalil 


+0 90888311 
=e 27266494 
4037753605 
~¢ 32360235 
+019416141 
-0 87372634 


J 15 16 


+026906065 
=010762426 
+0 20258684 
-« 24010293 


J 19 20 


+¢82206352 
-e41103176 


J 4 4 


+0 66666667 
= 053333333 
+017777778 
— 0 33862435 
+042328043 
= 037624925 
+025083285 
— 13030278 
+054292822 


52 
53 


50 
50 
50 
49 
48 
47 
46 
45 
43 
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Jos 5 
4081380208 
084771051 
4+037844220 
098552656 
4017109835 
021387294 
420253121 
~+ 15069286 
4290560614 


it, Oo aeo 


+¢20652090 


=042225482 © 


4038866635 
=021862484 
4085137559 
024629080 
+055415430 


J 213413 


4059373592 
~017918102 
4025234661 
022211758 
+013800688 


31717 


+017744069 
=e 71222722 
+0e13541689 


Seorenrl 


4012500000 
= 010416667 
+032552085 
= 054253476 
4$056514040 
040367172 
$2021024568 
= 83430820 
26072133 
—'e65838720 


wo) 6. 25 


+¢33908420 
= ¢30275375 
4¢11826319 
=—027375737 
+°042774587 
= 048607487 
+ 042194001 
= 028979396 
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BESSEL 


J 6 6 
+e10125000 
=e 13017857 
+e 73225455 
=—0 24408482 
+2054919081 
=e 89867592 
+¢11233449 
=e 11110006 


J 10 10 


+2026911444 
=—061162373 
+¢63710811 
=e 40840262 
+e18232261 
=~060774197 


J 14 14 


+¢77797246 
=e 25413766 
+e38914831 
-~¢37388759 


J 18 18 


4023443593 
=0 99943736 


Jae tre 


+0 16666667 
= 041666667 
4+ 041666667 
=~ 023148149 
4082671959 
020667989 
4038274054 
~054677221 
4062133204 
= 057530745 


J 7 6 


+043392857 
048816963 
4024408484 
0 73225447 
$014977932 
022466896 
$025923349 
023807155 


51 
51 
50 
50 
49 
48 
48 
47 


51 
51 
51 
51 
51 
50 


51 
52 
52 
52 


52 
52 


FUNCTIONS 


aaa Sera / 
+012765734 
—019547529 
+e13303180 
—054321318 
+015123548 
—¢ 30877244 
+048493109 
—2 60616389 


J 11 11 


+234900678 
=087978794 
+010235994 
-o 73723534 
+e37168946 
=e 14054509 


J 15 15 


+010219257 
=e 35927077 
+059438178 
061914771 


J 19 19 


+031020863 
~0 13998164 


3-453 


421093750 
— 094921875 
4017797852 
~019069127 
4013407980 
~°67039898 
4025139961 
_+ 73460930 
4017217406 


Jeet 


+o55850084 
=e 76018170 
+¢46561131 
—~017284056 
+°44110350 
~ 083131047 
4+¢e12123278 


52 
53 


J 8 8 
+e 16253968 
~« 28895943 
+e23116756 
—e 11208124 
+e 37360413 
—091964093 
+017516971 


J 1212 


+0 45444155 
—e 12584535 
+e°16180118 
—-0 12944094 
+e 72810527 


J 16 16 


+e 13453033 
—050646711 
+0 90038599 


J 20 20 
+¢41103176 


J 5 4 


+¢26666667 
-—¢17777778 
+¢50793651 
- 084656086 
+¢94062317 
-« 75249850 
+¢45605971 
021717129 
+¢83527418 


a FE) 


+e 12239856 
=«11558377 
+ 084060925 
= ¢ 37360414 
+¢11495511 
= 026275456 
+°46711922 


51 
51 
51 
51 
50 
49 
49 


52 
52 
52 


52: 


13 


14 


J10 9 


+093834405 
017274061 
+014574989 
=e 75677829 
+027365644 
-e 73887236 


J 14 13 


+°27566311 
—¢77645107 
+°10251581 
=¢ 84927309 


J 18-17 


4083791438 
=e 31862797 


SER 1 


+e¢10000000 
+e 10000000 
—«12500000 
+¢52083336 
-0 10850695 
+¢13563368 
-011302807 
+¢67278612 
-0 30035096 
+°10428854 
28969034 
+e 10000000 
—2 33333334 
+041666668 
—027777778 
4011574074 
—033068784% 
+268893298 
—0 10935444 
+e13669305 
+e11250000 
063281250 
+01423828) 
-017797852 


~ #0 14301845 


— 080447870 
+033519948 
—0e 10774270 
+¢27547848 
+013333333 
—e 10666667 
+¢35555556 
~e67724870 
+084656086 
-« 75249850 
+250166570 
—0 26069556 
+¢16275042 


THE ASTRONOMICAL JOURNAL 


BESSEL 
J 11 10 


%o122324675 
= 025484323 
+024504157 
=e 14585808 
+060774199 
=e 18991938 


J 15 14 


+036305380 
011118523 
4016023753 
=e 14540073 


J 19 18 


+e11104859 
=0 44974680 


SER 2 


+< 10000000 
+050000000 
+¢37500000 
+031249999 
+027343740 
+024609390 
+022558580 
+020947310 
+019638050 
+¢18547070 


+e20000000 


+e75000000 
+067708332 
+058051206 
+e51675620 
+¢47003650 
+043399290 


+040511700° 


+e 38131720 
+e 36126600 
+e25000000 
+o 33333332 
+265624993 
+054097229 
+249059590 
+e 45027340 
+041842020 
+0¢39243280 
+e37072300 


. te 32500000 


=~0 39062509 
+e76757818 
+.49057610 
+046928230 
+043258410 
+¢40430640 
+¢38079060 
+0 36090600 


FUNCTIONS 
J 12 11 


015996144 
037221797 
4040212835 
=e27031960 
+012776825 


J 16 15 


+e47902770 
015850181 
+°24765909 


J 20 19 
+014734910 


SER 3 


+010000000 
+015000000 
+018750000 
+021875000 
+024609370 
+027070318 
+°29326170 
+031420920 
+033384690 
+035239410 
++30000000 
+033750000 
+040781250 
+046578773 
+051746457 
+0 56446830 
+060786760 
+064837960 
+068651100 
+072263780 
++45000000 
+235000000 
+044062497 
+049250000 
$054167525 
+058669896 
+062854090 
+066778460 
+070485620 
+066250000 
+030703123 
048345704 
051649657 
+956499800 
+060815210 
+064858710 
+068666680 
+072275660 


52 


J 1312 


+e 20974226 
053933726 
+064720469 
—0 48540351 
+0 25697832 


J 17 16 


4063308388 
~022509649 
+¢37910988 


J 21 20 
+0 19572961 


SER 4 


+» 10000000 
+200000000 
+200000000 
+200000000 
+°40000000 
-»50000000 
~ «15000000 
+0 17000000 
+0 18000000 
- «20000000 
+000000000 
+200000000 
+¢00000000 
+0 60000000 
+240000000 
+¢90000000 
-+ 19000000 
~« 79000000 
- «70000000 
-. 20000000 
4000000000 
4000000000 
+000000000 
+000000000 
+e90000000 
-« 20000000 
+200000000 
~¢11000000 
+e 20000000 
+200000000 
+000000000 
+e 10000000 
-«40000000 
-e 23000000 
+0 20000000 
+e50000000 
+e 10000000 
+2 00000000 


52 


64, No. 126¢€ 


1959 February 


SER 1 


-—° 16954210 
+275688440 
-—219710531 
+ 234219670 
—042774588 
+040506242 
—» 30138572 
4020250000 
~e 26035714 


+014645091 


048816964 
+010983816 
~-017973518 
+022466898 
+¢25531468 
-039095058 
+e26606360 
—2 10864264 
+¢30247096 
061754488 
+¢96986218 
+032507936 
— 057791886 
4046233512 
—0 22416248 
+e74720826 
—e 18392819 
+041704180 
-@ 84450964 
+077733270 
—043724968 
+e17027512 
—e49258160 
+¢53822888 
-0 12232475 
+e12742162 
—081680524 
+0 36464522 
+069801356 
=e 17595759 
+¢20471988 
-0 14744707 
+¢74337892 
+°90888310 
—25169070 
+¢32360236 
-«258886188 
+°11874718 
—¢35836204 
+050469322 
— 044423516 
+2 15559449 
-—50827532 
+077829662 
+°20438514 
—0 71854154 
+e11887636 
+°26906066 
010129342 
+035488138 
—0 14244544 
+¢46887186 


SER 2 


+042916666 
—e 16125001 
+e11618056 
+e37987790 
+¢45986010 
+041592520 
+° 39153550 


+237008040. 


+057135417 
—e 36069878 
+0 21062479 
+090259300 
+°48952280 
+¢ 39588170 
+e 38038310 
+e 36017000 
+e 76437500 
—067758927 
4040555249 
-—063002220 
+263675410 
+e 35610990 
+e 37303540 
+e10258898 
=e 11706620 
+e 77613399 
—e22767056 
+0 10954314 
+024579980 
+e37882100 
+¢13799678 
=~ 19254719 
+e 14441040 
=057772094 
+0e22996709 
=e 73217600 
+e 18593040 
—e 30661132 
+026034198 
—0 12802456 
+051572729 
—0 94145700 
4+e25082627 
—0 47718759 
+0¢45590443 
—e 26277027 
+e11472865 
+e 33870215 
-—e 73008756 
+e 77852092 
—051202999 
+¢24696023 
+045772023 
=-e11023521 
+e13012491 
—095975672 
+0e61895163 
—0 16469844 
+021357317 
—-0 17446471 
+0e83741410 
=e 24396315 


E cos 


54 
56 
58 
60 
62 


SER 3 


4096249999 
+016124999 
+056187500 
+053213797 
+058828257 
+062882220 
+066805520 
+070505990 
+013851562 
016233724 
4072948352 
+052169362 
+061513130 
+064830410 
+068702990 
+072299150 
+219793750 
078107138 
+010889901 
4044018420 
+065763280 
+066448720 
4070583320 
4028130925 
018761248 
4018267087 
+018668390 
4074993310 
4066953210 
+072583890 
4039806027 
=037240579 
+032680674 
044792850 
+097757170 
+063918130 
+056126802 
067430390 
+059657066 
218732842 
4015371195 
+050707150 
¥0 78905516 
—e11560294 
+010837776 
248588103 
+028520885 
+011065171 
019110191 
+019383001 
—. 10803300 
4057905026 
+015483800 
030776420 
4$034012445 
022183860 
+021626650 
048597499 
+058557662 
043293489 
4030157369 
- 75562997 


—E COS 


54 
56 
58 
60 
62 


64 - 


66 
68 
55 
57 
59 
61 
63 
65 
67 
69 
56 
58 


54 
54 
54 
54 
54 
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SER 4 


=e 10000000 
=e 10000000 
- 30060000 
+e00000000 
+e 30009000 
+200000000 
+e 10000000 
- 90000000 
+e 10000000 
+e 20000000 
—e 10000000 
-e 30000000 
+e20000000 
-e70000000 
+240000000 
+e90000000 
=e 10000000 
—e 10000000 
~ «50000000 
+e 20000000 
+e 20000000 
+200000000 
—«60000000 
+200000000 
4000000000 
+e 10000000 
=e50000000 
4+e50000000 
=e90000000 
+e 13000000 
4+e00000000 
+e 19000000 
=e 80000000 
+211000000 
=e 10000000 
=e 10000000 
~e 10000000 
+200000000 


~e 18000000: 


—»29000000 
+* 10000000 
4080000000 
4000000000 
4000000000 
+0 40000000 


+e 26000000: 


+0 40000000 
+e 00000000 
-¢ 10000000 
-e 11000000 
-«e 60000000 
+e 70000000 
4000000000 
+e 70000000 
+e 70000000 
+e 10000000 
+e00000000 
+e 00000000 
—¢ 50000000 
+e 20000000 
+e00000000 
+e 20000000 
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SER 1 
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E cos 


+2062041726 52 69 69 


SER Gn 


-e 10000000 
+e 10000000 
-e11250000 
+013020833 
~053168406 
+010986329 
—e 13676396 
+¢11370086 
—067578963 
+ 230139385 
—0 10457823 
+e10000000 
= 213333333 
+037500001 
— 044444446 
+ 028935185 
—011904762 
+033757717 
—2069986842 
+e11072137 
013807379 
+¢11250000 
—017578125 
+¢77519531 
~e16018066 


+¢19228037 


—015106324 
+083799865 
~ 034597375 
+011049748 
+013333333 
—»24000000 
4014222223 
042328043 
_ 076190479 
092161071 
+080266507 
052772626 
+027146412 
+16276042 
033230252 
+024523054 
-095398971 
+e23132498 
~ 38497129 
+046825212 
~043520099 
+0 20250000 


SER 2 


+034513779 
+011334806 
— + 35880066 
+ 055036633 
+215347915 
= 052451533 
4020788493 
076281019 
+»28165427 
+9 38169287 


SER 6 


4210000000 
~ «87500000 
4088541667 
029405382 
~018488227 
—» 14234201 
011225149 
—291060063 
~° 75633187 
~264048879 
+2 10000000 
~° 11666667 
4033333334 
—e52777780 
~ 268865735 
—280348875 
—2 65289484 
054592860 
046309317 
—2 39833796 
4e11250000 
—e 16171875 
4071894533 
~e 15798341 
4012385125 
~° 66235661 
~e43371171 


37899070 . 


032721582 
$0 13333334 
022666668 
$013444445 
040820106 
4069080691 
012570914 
024465817 
~» 28727076 
—2 24746399 
4016276042 
—. 31873915 
4023457361 
=091985422 
+#22025908 
=~ 40049944 
4020797916 
=e 26228653 
4020250000 
244839286 


E cos 


68 
65 
67 
69 
66 
68 
67 
69 
68 
69 


SER 3 


+099044170 
+041992578 
-»11603955 
+016486217 
+ 058397619 
17638207 
4081118328 
~ 026580903 
+011256236 
4015604857 


SERT 


4075000000 
+012500016 
+ 046874700 
+023437800 
413671600 
4087902000 


+060420000 


+¢43660000 
+032740000 
-»20000000 
$018333334 
~026041720 
+065190800 
034078200 
4$221353200 
+014333400 
4210156000 
+e 75020000 


+e57000000. 


+0 10000000 
«40000000 
4034791664 
—»69861080 
+011822860 
+020583200 
4016926400 
+011732000 
+©85600000 
4-0 64480000 
+¢20000000 
-«67500000 
+064968750 
022296874 
+043450360 
-024767400 
+16508600 
+095360000 
+¢72560000 
+056220000 
+032500000 
=»10791666 
+011668055 
— 054193448 
+014190576 
-021096220 
+035884000 
+260660000 
+ 263420000 


—E cos 


SER 4 


+e12000000 
+e 00000000 
—e20900000 
—« 30000000 
+e 10000000 
+«20000000 
+°00000000 
+e 10000000 
+2e00000000 
+¢00000000 


SER 8 


-« 20000000 
+e 16666666 
-e 78124800 
+e52429900 
+0 34798300 
+e 26760800 
+e21260000 
+e 17368000 
+e 14511000 
+e 12372000 
+e 10000000 
-e 40000000 
+¢33958335 
—e 73194460 
+e 10399290 
+e 15737900 
+e 16447000 
+e 13301000 
+e 11082000 
+e 93990000 
+e 20000000 
-—2 67499999 
+064406246 
—e22507810 
+¢ 42257540 
-—e¢ 30232000 
+e 14427500 
+e 92910000 
+e 80290000 
+e69060000 
+e 32500000 
—e 10791667 


+e 11623611 | 


054320439 
4014091379 
-» 21597390 
+2 33404000 
4051320000 
+» 63000000 
+049166667 
-e 16776042 
+0 20219556 
011648544 
+0 38964414 
—e 83610490 
+e 13508600 
-« 89020000 
+¢61890000 
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SER 5 


046285714 
+040680805 
-»219526787 
+259800780 
—012781168 
4020220208 
024688899 
$025531468 
064626526 
- 4065701418 
037470624 
+013888974 
036422545 
+071453110 
4032507936 
290299822 
4010402540 
~ 068649760 
+029888331 
093113645 
4+021896213 
4041704180 
012615514 
4016218423 
012145824 
4060752480 
-021953328 
4053822888 
017614764 
4024974637 
020910214 
4011814505 
048619361 
4069801356 
~e 24575895 
4038067747 
~» 35216695 
+022178496 
+090888310 
~e 34257901 
4057529306 
058248424 
4040450293 
4011874718 
047710922 
+086305526 
094892838 
4015559449 
~966386981 
4012865719 
-015260718 
4020438514 
092292668 
+019073051 
+026906066 
012819949 
4028137062 
+0 35488138 
~0 17793358 
4046887186 
024677466 
4062041726 


THE ASTRONOMICAL JOURNAL 


SER 6 


+0 39234375 
18893974 
+057721540 
012635292 
$017737760 
—041792721 
4025531467 
063030809 
4063751100 
~+ 36394069 
+0 13487760 


we 35662225 


4068160292 
032507937 
~e 28493827 
+210140670 
~0 66905331 
$029134749 
—0 91083229 
$021284441 
+#041704179 
=. 12406993 


+0e15867730 ! 


—e11871992 
+059387419 
—e 21497105 
+053822889 
—017370115 
+024505724 
020489918 
4011575367 
—¢e 47678783 
+069801355 
=0 24285055 
+037441326 
—e 34582334 
+e21771300 
+¢90888311 
—e 33908333 
+0 56692832 
=—0 57303333 
+e 39773158 
4¢¢11874718 
—0 47286825 
+e85188735 
—0 93499337 
+¢ 15559449 
~e 65868335 
4e12716608 
—0 15056977 
4020438516 
—0 91653966 
+e 18873926 
+0e26906065 
—-° 12740813 
4+e27871078 
+e 35488138 
—0 17694779 


4046887183 


0 24554079 
4062041727 
$082206352 


E SIN 


59 
61 
63 
65 
67 
69 
56 
58 
60 
62 
64 
66 
68 
57 
59 
61 
63 
65 
67 
69 
586 
60 
62 
64 
66 
68 
59 
61 
63 
65 
67 
69 
60 
62 
64 
66 
68 
61 
63 
65 
67 
69 
62 
64 
66 
68 
63 
65 
67 


56 
56 
56 
56 
56 
56 
57 
57 
57 
57 
57 
57 
57 
58 


SER 7 


+049166666 
016776041 
4020258308 
011642489 
4039052036 
083182020 
+013755200 
=e 76880000 
4066200000 
+071875000 
=025612500 
034229798 
023223214 
#095126152 
=e25836020 
4051427200 
— 68636000 
4010305208 
~0 38598220 
056586388 
044023758 
021264768 
~069866320 
4016795020 
229948200 
4014596181 
057588218 
4091890152 
080368116 
+044634930 
~017174986 
+048752460 
4020504408 
085232366 
+014702531 
014248067 
4089291772 
239335590 
4012903614 
+028636098 
012530897 
4023232262 
024671248 
4017197669 
085274522 
+039819888 
—0 18319479 
4$036320864 
041896510 
+032119424 
017690945 
+055188722 
=0 26652692 
056261748 
-069992770 
4058478586 
+076291452 
~2 38613384 
+086450830 
011530476 
+010420561 
4010524502 


—E cos 


SER 8 


+e 71875001 
— 025612500 
+0 34193639 
—¢ 23223214 
+095040256 
—e 25870070 
+051191000 
-e69915000 
+e 10305209 
—0 38598221 
+056550929 
=044017557 
+4e21255215 
~0 69889440 
+e 16772370 
-« 30072100 
+e 14596181 
-0 57588220 
+2091854033 
—« 80354980 
+044622953 
—e 17175743 
+e 48729240 
+e 20504409 
—e 85232367 
+o 14698740 
—0 14245935 
+¢89275481 
=—e 39333898 
+e 12900919 
+e 28636099 
—o 12530897 
+e23228186 
-2 24668112 
+e 17195348 
=e 85268915 
+e 39819888 
—2 18319479 
+e 36316390 
=e 41892116 
+e 32116050 
—6 17689754 
+055188722 
—e 26652691 
+e 56256756 
—069986778 
+o 58473658 
+e 76291453 
—¢ 38613384 
+e 86445177 
—-e 11529672 
+e 10419843 
+2 10524503 
=e 55734403 
+e 13189174 
=—0 18764986 
4-0 14494207 
=080182290 
oe 19996342 
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SER 5 E cos SER 6 E SIN SER 7 E cos SER 8 E SIN 
+e82206352 52 69 70 —055734402 53 64 64 +¢30145968 54 68 66 
4013189822 54 66 64 +227381731 53 65 67 

—018766054 54 68 64 —016454258 54 67 67 

+014494207 53 63 65 +045216714 54 69 67 

—e80182288 53 65 65 +¢37575981 53 66 68 

+019997093 54 67 65 -«23483332 54 68 68 

030219332 54 69 65 +4051521884 53 67 69 

+e19933403 53 64 66 —0«33442012 54 69 69 

-011501493 54 66 66 +2e70591985 53 68 70 

5 +296659047 53 69 71 

SER 9 — cos SER 10 E cos SER ll E cos SER 12 E SIN 
+e10000000 51 50 50 +.10000000 51 50 50 50000000 50 51 50 +210000000 51 50 51 
+e50000000' 50 52 50 +.215000000 51 52 50 42e10000000 51 50 51 -—e12500000 50 52 51 
-e10000000 51 51 51. =e20000000 51 51 51 -—0e37500000 50 52 51 +0¢52083336 48 54 51 
+037500000 50 53 51 +4¢25000000 50 53 51 +40¢26041667 49 54 51 -010850695 47 56 51 
026041667 49 55 51 ~010416667 49 55 51 275954865 47 56 51 +4023563368 45 58 51 
+2075954865 47 57 51 +42e21701390 47 57 51 4012207032 46 58 51 -—e11302807 43 60 51 
—e12207032 46 59 51 —027126736 45 59 51 012433089 44 60 51 +4067278612 40 62 51 
+012433089 44 61 51 4022605614 43 61 51 +40¢87462205 41 62 51 —«e30035096 38 64 51 
—087462205 41 63 51 —013455722 41 63 51 —e45052645 39 64 51 +0¢20428854 36 66 51 
+045052645 39 65 51 +4+0e60070192 38 65 51 4017729049 37 66 51 —«28969034 33 68 51 
—017729049 37 67 51 1020857707 36 67 51 055041166 34 68 51 +«50000000 50 51 52 
4055041166 34 69 51 4057938068 33 69 51 450000000 50 51 52 -+16666667 50 53 52 
-—050000000 50 52 52 250000000 50 52 52 33333334 50 53 52 +4e20833334 49 55 52 
+033333334 50 54 52 +4¢16666667 50 54 52 4462500002 49 55 52 —-13888889 48 57 52 
—062500002 49 56 52 020833334 49 56 52 ~655555557 48 57 52 4257870371 46 59 52 
+055555557 48 58 52 +0613888889 48 58 52 4028935186 47 59 52 -216534392 45 61 52 
— 028935186 47 60 52 ~057870371 46 60 52 099206350 45 61 52 +034446649 43 63 52 
+°99206350 45 62 52 4016534392 45 62 52 +4024112655 44 63 52 -e54677221 41 65 52 
-024112655 44 64 52 634446649 43 64 52 =043741776 42 65 52 +068346524 39 67 52 
+°43741776 42 66 52 +4+0e54677221 41 66 52 068346510 39 67 52 -—269036893 37 69 52 
+0¢68346510 39 68 52 -—068346524 39 68 52 069036900 38 69 52 4237500000 50 52 53 
-°37500000 50 53 53 -—«e25000000 50.53 53 +e37500000 50 52 53 -—.21093750 50 54 53 
+¢35156250 50 55 53 +2¢14062500 50 55 53. ~e35156250 50 54 53 +¢47460937 49 56 53 
-e11074219 50 57 53 —031640625 49 57 53 4011074219 50 56 53 —.59326173 48 58 53 
+017797852 49 59 53 +¢39550782 48 59 53 -e17797852 49 58 53 +.647672818 47 60 53 
-°17480033 48 61 53 —~e31781878 47 61 53 +4€17480033 48 60 53 226815956 46 62 53 
+°11620250 47 63 53 +617877304 46 63 53 ~0e11620250 47 62 53 4011173316 45 64 53 
—0©55866582 45 65 53 =—274488773 44 65 53 +255866582 45 64 53 235914232 43 66 53 
4020351398 44 67 53 +¢23942821 43 67 53 ~e20351398 44 66 53 4091826159 4il 68 53 
058156569 42 69 53 =—261217439 41 69 53 +4e58156569 42 68 53 4033333334 50 53 54 
- 033333335 50 54 54 -—016666667 50 54 54 +033333335 50 53 54 -e26666667 50 55 54 
+¢40000002 50 56 54 +613333333 50 56 54 —e40000002 50 55 54 +¢88888890 49 57 54 
-el7T777779 50 58 54 —044444445 49 58 54 4617777779 50 57 54 016931218 49 59 54 
_ +042328043 49 60 54 +084656088 48 60 54 —042328043 49 59 54 +621164022 48 61 54 
263492067 48 62 54 =010582011 48 62 54 +063492067 48 61 54 -018812463 47 63 54 
+ 065843622 47 64 54 +40694062313 46 64 54 —265843622 47 63 54 +012541643 46 65 54 
050166568 46 66 54 =062708213 45 66 54 +450166568 46 65 54 -—e65151390 44 67 54 
+029318124 45 68 54 +4032575695 44 68 54 ~.29318124 45 67 54 +627146411 43 69 54 
—032552084 50.55 55 =—.13020833 50 55 55 4013573206 44 69 54 +632552083 50 54 55 
047471788 50 57 55  +013563368 50 57 55 +4032552084 50 54 55 -e33908420 50 56 55 
027247837 50 59 55 060550752 49 59 55 247471788 50 56 55 +615137688 50 58 55 
+086726340 49 61 55 +4015768425 49 61 55 +027247837 50 58 55 -—e39421062 49 60 55 
017794229 49 63 55 =—027375736 48 63 55 086726340 49 60 55 +668439340 48 62 55 
+025664752 48 65 55 +4¢34219670 47 65 55 +42017794229 49 62 55 —.285549176 47 64 55 
027544243 47 67 55 —032404994% 46 67 55 025664752 48 64 55 +¢81012484 46 66 55 
+#022905315 46 69 55 +024110858 45 69 55 +4¢27544243 47 66 55 -—060277144 45 68 55 
-¢33750001 50 56 56 011250000 50 56 56 022905315 46 68 55 +233750001 50 55 56 
4057857144 50 58 56 4614464285 50 58 56 4633750001 50 55 56 —24339285# 50 57 56 


1959 February 


| SER 9 


} 


| 


= 040680805 


+016272322 
042714842 
+079882307 
011233449 
~ 036473524 
+071807249 
~ 059728563 
+ 028823555 
292593154 
4021425026 
- 037606902 
~ 240634921 
+ 090299622 
- 086687832 
049035542 
— 218680206 
+051729798 
046337976 
+011468650 
=012475711 
+080972158 
- 035736751 
+ 011554383 
= 053822889 
+e14678970 
-.17839025 
+013068884 
= 065636136 
~ 063455777 
+0 18904534 
025378499 
4020715703 
011672893 
~275740259 
4024469930 
-035955814 
+032360234 
091343985 
4031807282 


=050767956 


+2049943600 
—e 11113892 
+041491864 
071476217 
—0 13625677 
+054289806 
=. 10038448 
-» 16816291 
4071221937 
—e 20875375 
+093649249 
-0 26048435 
032653541 
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SER 10 


081361616 
4027120535 
061021201 
4099852879 
=212481610 
010421007 
+015957167 
=010859739 
+044343933 
=012345753 
4025205913 
=. 39586211 
=010158730 
+0 18059964 
014447973 
+070050775 
023350258 
+057477558 
=. 10297328 
+020852090 
=019193400 
+010796288 
042043239 
4012162509 
010764578 
024464949 
025484324 
+e16336105 
072929044 
011537414 
+029083899 
=o 33837996 
024371416 
012287255 
012623376 
+0 34957042 
044944773 
4035955817 
=014052921 
042409709 
=059727008 
+052572208 
= 15876989 
4051864828 
=079418021 
=.18167568 
+063870360 
—. 10566787 
=021020364 
+079135486 
= 024559264 
+098578159 
= 028942707 
= 034372146 


SER 11 


057857144 
+040680805 
016272322 
$042714842 
—079882307 
011233449 
~012344451 
036473524 
071807249 
$059728563 


—028823555 - 


+092593154 
021425026 
+037606902 
+040634921 
=—090299822 
+e86687832 
—049035542 
+¢18680206 
-—051729798 
+010948106 
+ 046337976 
—011468650 
+012475711 
=-«080972158 
+035736751 
— 011554383 
+053822889 
-—014678970 
+017839025 
—e 13068884 
+065636136 
—024309681 
+063455777 
—0 18904534 
+025378499 
-—020715703 


- +011672893 


+075740259 


3024469930 


4035955814 
032360234 
$020225146 
+091343985 
031807282 
+050767956 
049943600 
+#011113892 
041491864 
+071476217 
=0 76303589 
4013625677 
= 054289806 
+010038448 
+e16816291 
071221937 
+014068531 
4020875375 
093649249 
+026048435 
~e 12338734 
$032653541 
=+ 16441270 


SER 12 


+0 24408485 
081361608 
+0 18306361 
=0 29955865 
$037444831 
~ 037033354 
$0 36473525 
055850082 
+0 38009085 
015520376 
4$043210136 
~0 88220695 
$0 13855174 
+040634920 
072239858 
4057791890 
=0 28020310 
+#093401033 
022991023 
4043792428 
+046337977 
093834404 
+086370299 
0 48583297 
+0 18919457 
054731288 
4053822888 
012232475 
+012742162 
081680524 
+0 36464522 
=012154839 
+063455778 
015996144 
+018610898 
013404279 
+067579901 
+0 75740260 
=0 20974225 
026966864 
21573490 
+#012135088 
+091343985 
=. 27566310 
+0 38822554 
034171935 
#011113892 
—+ 36305379 
+055592615 
053412512 
+013625676 
047902768 
+#079250902 
016816291 
063308389 
$011254825 
020875376 
083791438 
$0 26048436 
011104859 
4032653541 
40 16441270 
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The general formula was then arranged in the following form: 


Ix(ie) = (2) az z Gig (eas a ( 


aes ae 


x 


x ( 


The calculations begin at the left end of the 
second formula, and each X signifies a separate 
multiplication operation. It is apparent that the 
punched results of each pair of operations are 
used immediately as input for the next operation 
until the final term on the right is reached. Each 
time the start was made at 2m + k = 20 or 21. 
Finally, we applied the check equation: 


ees fee = < [Teqlee) age). 


The potentialities of this new technique are 
illustrated by the results given in the adjoining 
tables. These are to be read as follows: under 
the Series No. appears the numerical coefficient, 
followed by its exponent of 10 (on the base of 
50 for zero). The index under £ is the power of 


—7e — 7? I I te \! 
Pe )+axml() 


the eccentricity, and the index under COS o 
SIN is the coefficient of the mean anomaly of : 
cosine or sine term. Again both indices are give: 
to the base 50 for zero. Each Bessel Function 
J;,(ie), carries a heading J k 2. Series No. 1 i 
(a/r), No. 2 is (a/r)?, No. 3 is (a/r)’, No. 4 i 
(a/r)?(r/a)? = 1 (for a check), No. 5 is cosf 
No. 6 is sin f, No. 7 is cos2f, No. 8 is sin 2. 
(f = true anomaly), No. 9 is (7r/a), No. 10 i 
(r/a)?, No. 11 is cos X, No. 12 is sinX (X = 
eccentric anomaly). Series No. 7 is slightly cur 
tailed for typographical reasons. All the coeffi 
cients were derived from the Bessel Functions 


eg., @/r =I+ 3 2J;(ie) cos 7(1), or the serie 


t=1 
were developed from each other, e.g., sin 2f = 
2 sin f cos f. 


A GRAPHICAL METHOD FOR THE DETERMINATION OF THE 
EPHEMERIS FOR VISUAL DOUBLE STARS 


_ By D. BELORIZKY 
Observatoire de Marseille, Marseilles, France 


Received May 1, 1958 


Abstract. A simple representation is given for the eccentric anomaly in the apparent motion of double stars. A graphica 


method is deduced for the determination of the ephemeris. 


__ 1. The calculation of the ephemeris for visual 

double stars makes use of the Thiele-Innes for- 
_mulae. Let E be the eccentric anomaly of the 
real orbit, e its eccentricity; let 


X=cossE—e, Y=Vi-—ésink 
one has 
x =pcosd=AX+ FY; 


y = psind = BX +GY, (1) 


where p is the apparent distance, and @ the posi- 
tion angle of the companion star in a couple. 


This equation, taken together with that o: 
Kepler 


E—esinE = p(t — T), (2: 


where uw = 27/P, P the period, T the time o 
passage at the periastron, yields all the quantities 
necessary for the calculation of the ephemeris. 

Equations (1) are due to Innes (1926). Thiele 
(1883) had given them much earlier unde 
another form, viz. 


x = a,cos p(cosE —e)+bicosgsinE (3) 
y = asin p(cos E — e) + bi sin gsin £, 


1959 February 


90° 
y 


Figure 1. The apparent orbit. 


where a; is the half-diameter of the apparent 
orbit passing through the principal star I, b1 the 
conjugate half-diameter, p the position angle of 
a, and g the position angle of dy. 

One sees at a glance that 


A=aicosp, B=aisinp; 
pa Pos g _ Sising (4) 
ee Ape 


Without disparaging Thiele’s and Innes’ con- 
tributions in any way, it is in order to note that 
Lagrange had used, in a note read before the 
Institut de France on August 22, 1808, the ec- 
centric anomaly instead of the true anomaly, in 
order to get the projection of an elliptic motion 
on any plane (Lagrange 1815). Moreover this 
mode of projection of an elliptic motion is always 
used in the theory of. variation of arbitrary 
constants (Tisserand 1889). 

2. The author has obtained formulae similar 
to (3), making use of the apparent orbit only 
(Belorizky 1948). 

Take the conjugate diameters OP and OR (see 
figure) of respective lengths ai, bi, as x1 and 41 
axes : the orbit can then be given by the equations 


%1 = a,cosu, yi = bisin u, (5) 


where u is a parameter. 

Let ¢: be distance IO, I being the principal 
star: it is well known that c;/a; = e. The law of 
areas yields the relation 


d dx 
(x1 — en - Wi ae = 01b1; 
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taking into account (5) and integrating, one gets 
(6) 


On the other hand, if the origin of co-ordinates 
be shifted to I and rectangular axes Ix, Iy be 
taken (Ix pointing North) one has formulae (3), 
with the parameter u replacing eccentric anomaly 
E. Therefore u = E. One thus finds in the ap- 
parent orbit a remarkable and very simple 
representation of the eccentric anomaly of the 
real orbit. 

3. Formulae (5) and (6) yield a graphical 
method for the determination of the ephemeris of 
double stars. Let a half-circle be drawn with OQ 
on diameter; draw chord OA. Angle QOA = u; 
take OM = OA along OP; let MC be parallel to 
OR: then CI = p, and angle NIC = @; equation 
(6) gives ¢. Varying u by steps of 10° for instance, 
one has the ephemeris. When u > 90°, let 
u = 90° + wu and x, be taken equal to —QA, 
angle QOA being equal to 1. 

Conversely, if » and 6 be measured at any 


u—esinu = p(t — 7). 


time, the worth of the measurement can be 


checked without having to draw the apparent 
orbit: CM being taken parallel to OR and OA 
equal to OM, one has uw. Equation (6) gives ¢, 
to be checked against the time of observation. 
On the other hand, if C be on the orbit, 
(CM/OR)? + (OM/OP)? should be equal to 1. 

Note finally that the geometrical elements 
needed for the above graphical construction are 
a1, bi, p, gq; they appear to be easier handled than 
A, B, F, G, for a; and 0, are intrinsic elements of 
the ellipsis, p and g conventional elements de- 
pending on the choice of axes; while (4) shows 
that A, B, F, G are functions both of a, 6; and 
p, gq. These elements are thus not well adapted 
for the calculation of the elements of the orbit, 
and particularly for that of the half-major axis. 
Moreover, one has to turn back to di, 01, p, and g 
for the construction of the apparent orbit. These 
last quantities are more convenient for the de- 
termination of the elements of the real orbit 
(Belorizky 1957), and appear quite naturally in 
that of the apparent orbit. 
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ON THE. COMPUTATION OF STELLAR SPACE DENSITY FUNCTIONS 


By HAROLD K. CROWDER 
Case Institute of Technology, Cleveland, Ohio 


Received September 17, 1958 


Abstract. A program has been developed for computation of space density functions on the IBM 650 magnetic drum 
calculator. Output, based on an analytic solution of the fundamental equation of stellar statistics, is compared with data 
previously obtained by the (m, log x) method. In regions free of conspicuous irregularities in density the results of the 
two methods agree sufficiently well to warrant extensive use of machine computation, taking advantage of large savings 
in time. For stars of spectral type F8 to K3, densities obtained by any method can only be approximations because of 
uncertainty as to the true percentages of giants and dwarfs in these categories. 


McCuskey (1956) has published a summary of 
a sequence of papers setting forth the results of 
analysis of star counts made at the Warner and 
Swasey Observatory of Case Institute of Tech- 
nology. Included in the summary are tables of 
stellar space density functions for stars in selected 
regions of the Milky Way according to spectral 
types. Centers of the regions are listed in the first 
paper of the sequence (McCuskey and Seyfert 
1947). Computation of the space densities was 
time-consuming; and, upon the establishment at 
Case of a computing center equipped with high- 
speed electronic computing machines, it appeared 
desirable to construct a program for their ma- 
chine computation. 

The first step in the determination of the space 
densities is the solution of an integral equation 
to obtain a fictitious density function, which 
must then be corrected for interstellar absorp- 
tion. The equation is 


A'(m) = A(m)dm = ile ¢(M)dM - A(ro) dV, 


where 


m = apparent magnitude, . 
ro = apparent distance from the sun in 
parsecs, 
M = absolute magnitude 
=m-+5— 5lognro, 
dV = element of volume in thousands of 
cubic parsecs, 
A(ro) = fictitious density function = num- 
ber of stars per thousand cubic 
parsecs in dV, 
¢(M)dM = fractional part of the stars in dV 
having absolute magnitude in the 
range M — 3dM to M + 3dM, 
A'(m) = A(m)dm = number of stars per 
100 square degrees observed in the 
range m — dm to m + 4dm. 
It is assumed that ¢(M) is Gaussian, 


I nee (M—M,) 


Del ca a ’ 


where 


M, = mean absolute photographic magnitude 
per unit volume of space, 
o, = dispersion. 


The fictitious density function A(ro) was de- 
termined by the (m, log 7) method (Bok 1937), 
a trial-and-error method of numerical integration 
employing equal increments in log 79. Bok (1937) 
had noted that this method should be more sen- 
sitive to local clustering than analytic methods 
which had been proposed. However, it suffers 
the serious drawback that, as log 7o increases, the 
element of volume becomes so large that not only 
is local clustering not apparent, but the densities 
are averaged over such large volumes that they 
become almost meaningless. 

It was shown by Schwarzschild (1912) that if 
A(ro) has the form 


A(ro) = ebth log roti lox* ro, 


the integral equation can be‘directly integrated 
to yield a parabolic curve, 


log A’(m) =a + bm + cm’, 


where a, b, c are functions of M,, o,, h, k, 1. These 
equations can be solved for h, k, 1, and hence, if 
the log A’(m) curve is approximated by a parab- 
ola, the function A(ro) can be determined. The 
success of this method obviously depends on how 
closely the log A’(m) curve can be approximated. 
If the approximation is good the method pro- 
vides a point function for A(7o) not requiring 
averaging over large volumes of space. 

The correction for interstellar absorption was 
made by the method of Seeliger (Bok 1937), 


3h”) ’ 
Dir) = acres] + LO], 
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where 


r = true distance from the sun in parsecs, 
f(r) = total interstellar absorption in magni- 
tudes to distance 7, 
f' (r) = derivative of f(r) with respect to 7, 
mod = loge, 
D(r) = true density function in stars per thou- 
sand cubic parsecs. 


A program for the determination of A(7o) by 
the Schwarzschild method followed by the Seeli- 
ger correction has been prepared for the IBM 650 
Magnetic Drum Calculator using the Bell Float- 
ing Decimal Interpretive System (Wolontis 
1956). Input consists of a standard Bell System 
deck, a program deck, and a deck of 14 cards 
containing f(r) for r = 100 (100) 3000 parsecs 
(r ranging from 100 to 3000 parsecs in increments 
of 100 parsecs), the stellar data (1, and o,), the 
values of log A’(m) (base 10) using dm = 3, and 
certain information for finding the parameters of 
the approximating parabola. This information 
includes minimum and maximum values of m, 
which may be set at will as long as at least six 
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Figure 1. Log A'(m) and ZO curves for 
Bs stars in LF1 
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Figure 2. Tog A'(m) and D(r) curves for 
8-Ao stars in LF1. 
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values of m are used. Output includes the re- 
siduals between the log A’(m) curve and the 
computed approximate parabola, and the values 
of D(r) for r = 25 (25) 200 (100) 2900 parsecs, 
the computation being discontinued when the 
distance corresponding to the faintest observed 
magnitude (12.5) is attained. Provision is made 
for continuous computation of D(r) for different 
spectral types of stars and for different regions 
of the sky. For a fixed star type and region, M, 
and o, may be varied in order to study the effect 
on D(r) of these parameters. 

In order to test the analytic (Schwarzschild) 
method as compared with the (m, log 7) method, 
D(r) has been computed for all star types listed 
in the McCuskey summary for regions LF1, LF2, 
and LF4. Typical results, those for star types Bs, 
B8—Ao, A2—A5 and Fo-F5 in region LF1, are 
shown in Figures I to 4, respectively. The broken 
ascending curve is the observed log A’(m) curve, 
and the solid ascending curve is the approxi- 
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Figure 3. Log A’(m) and D(r) curves for 
A2-As5 stars in LF1. 
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Figure 4. Log A’(m) and D(r) curves for 
Fo-F5 stars in LF1. 
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mating parabola. The broken descending curve 
is the D(r) curve as computed by the (m, log 7) 
method, and the solid descending curve is that 
computed by the analytic method. 

In each of the illustrated cases the parabola 
provides a rather good approximation to the ob- 
served log A’(m) curve, although a few of the 
residuals (up to .096) are somewhat larger than 
the value (.05) adopted as an upper limit in the 
(m, log x) computations. In spite of this latter 
fact the correspondence between the D(r) curves 
computed by the two methods is entirely satis- 

‘factory. There is some discrepancy at relatively 

short distances from the sun, particularly evident 
in the B5 and B8—Ao stars. At these short dis- 
tances any analysis is questionable since the num- 
ber of stars involved is so small that the star 
counts used are not for specific regions but for 
the sky as a whole. In the case of the A2—A5 
stars the parabolic approximation is extremely 
good and the two D(r) curves differ only slightly 
from each other. Even in the case of the Fo—-F5 
stars, where the observed log A’(m) curve has a 
reversal of curvature, the two methods yield 
essentially the same D(r) curve. 

In an attempt to get a closer approximation to 
the observed log A’(m) curve it was felt that 
two parabolas might be used, one for the brighter 
stars and one for the fainter ones. Each parabola 
would yield a D(r) curve, one for stars close to 
the sun and one for those farther away. Because 
of the variation in absolute magnitude of the 
stars of a given spectral type there would be a 
region of uncertainty in trying to join the two 
curves, but by requiring an overlap of three 
points in fitting the two parabolas and then 
associating with the middle common apparent 
magnitude the distance to the star of most fre- 
quent absolute magnitude, a joining point could 
be determined. A slight smoothing in this region 
should then provide a better D(r) function than 
one developed through use of a single parabola. 
Many cases were computed on these two bases. 
The results indicated that the added precision in 
curve fit did not increase the credability of the 
resulting D(r) function. In fact, in a majority 
of cases the two curves were almost indistin- 
guishable. The gG8—K3 stars in LF1 provide an 
example. The residuals in the one-parabola and 
two-parabola fits are shown in Table I and the 
corresponding D(r) curves in Table II. 

In the cases of the F8 to K3 stars the analyses 
by the analytic and (m, log 7) methods do not 
agree very well. Here the observed log A’(m) 
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TABLE 1. [OBSERVED LOG A’(m) | — [PARABOLIC APPROXI- 
MATION | FOR gG8—K3 STARS IN LFI 


™m One parabola Two parabolas 
4.0 20 + .041 
4.5 — .078 — 028 
5.0 ++ .002 — .030 
5.5 +.051 —.037, 
6.0 +.129 — .013 
6.5 +.168 +.049 
7.0 +.116 +.021 +.041 
7.5 + .034 —.01I2 +.012 
8.0 — .048 —.018 —.028 
8.5 = 101 — .049 
9.0 = .123 — .052 
9.5 +086 — .007 
10.0 — .039 + .038 
10.5 =1,003 +.061 
II.0 — .016 + .022 
11.5 + .000 + .003 
12.0 + .027 — .019 
12.5 + .083 = 2021 


TABLE II. D(y) FUNCTIONS FOR gG8-K3 STARS IN LFI 


r One parabola Two parabolas 
50 -99 I.1I 
100 -44 -63 
200 -21 -19 
300 -14 ee 
400 II 09 
500 .09 07 

.08 

700 -07 06 
800 -06 06 
goo -06 06 
1000 .06 06 


curves are frequently highly irregular. Our stud- 
ies indicate that this irregularity may result 
largely from the uncertainty which exists in the 
relative percentage of giants and dwarfs among 
these spectral types (McCuskey 1956). In the 
light of this uncertainty it would appear that 
analyses of the F8 to. K3 stars cannot be as 
reliable as those of other types. Even in these 
cases the orders of magnitude of the results ob- 
tained by the two methods are consistent. 

The time element becomes paramount in the 
comparison of two methods which yield such 
similar results. It takes an experienced computer 
about six hours to determine by the (m, log x) 
method the D(r) function for a single star type. 
Using the electronic computer program devel- 
oped for the analytic method, D(r) functions 
were computed for I1 star types in LF2 in 
approximately 1} hours, the first hour being de- 
voted to key-punching the input data. The total 
machine time required was 24 minutes. A sec- 
ondary advantage of the analytic method is that 
the original values of log A’(m) may be taken 
directly from the star counts without preliminary 
smoothing. 
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In conclusion it may be stated that for the 
large majority of cases the machine computation 
for determining space density functions, using a 
single parabola to approximate the observed 
log A’(m) curve, yields results of sufficient v:Jid- 
ity to warrant extensive use of this method in 
‘regions characterized by relatively smooth log 
A'(m) curves. 

The author will be glad to furnish on request 
a program deck and full instructions for its use 
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to anyone having available an IBM 650 com- 
puter. 
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NEW PHOTOGRAPHIC METEOR SHOWERS 


By RICHARD E. McCROSKY ann ANNETTE. POSEN 
Harvard College Observatory, Cambridge, Massachusetts 


Recewed July 10, 1958 


Abstract. Seven new photographic meteor streams have been observed among 2538 meteors photographed with the 
Baker Super-Schmidt cameras. Approximate orbital elements were derived for each meteor. The average orbital elements 
‘of the new shower members were compared with known comet orbits, and associations were found between Comets 
1739, 18191V and 19131 and three new streams. Only two of the streams are readily identified in catalogues of visual 


meteors. 


The existence of a large number of weak meteor 
showers has been suggested by several observers 
(Denning 1899, Hoffmeister 1948, McIntosh 
1935) on the basis of their numerous visual ob- 
servations. Until recently the quantity of photo- 
graphic meteor data has, in spite of its high 
quality, been far too small either to confirm or 
to contradict this hypothesis. In the period from 
February 1952, through July 1954, 3623 meteors 
were photographed from two stations by the 
Harvard Meteor Project with the Baker Super- 
Schmidt cameras. A single visual observer prob- 
ably could not adequately observe and record 
meteors at more than twice the acquisition rate 
of the cameras. For the slowest meteors the cam- 
eras photograph to the visual limit; thus, quan- 
titatively, the photographic method is almost the 
equivalent of the visual. If most of the more 
prominent visual showers could be detected in 
the photographic material, this fact would con- 
firm the accuracy of the visual data. 

Of the meteors photographed during the above 
mentioned period, 357 have been accurately re- 
duced by Jacchia (unpublished) ; 2181 were re- 
duced by an approximate method (McCrosky 
1957) to yield velocities and radiants with mean 
errors of less than 3 percent and 3 degrees, re- 
spectively. These radiants, although far from 
representing the ultimate accuracy available with 


the photographic technique, are probably more 
accurate than visual radiants, and the addition 
of velocity as a second parameter in determining 
shower members improves the quality of photo- 
graphic data immensely. 

Orbital elements were computed for each me- 
teor with an IBM 650 program prepared earlier 
to determine elements of meteors reduced by 
accurate techniques (Whipple and Jacchia 1957). 
For expediency and economy, we used this pro- 
gram without change even though it contained 
refinements made superfluous by our compara- 
tively large errors. 

The radiants, corrected for zenith attraction 
and diurnal effect, were plotted on equal-area 
projections, one for each month of the year. The 
familiar major showers were readily identified 
and were analyzed to indicate the order of pre- 
cision of the orbital elements. The shower mem- 
bers were chosen by eye from the radiant plots. 
The radiants of the accepted shower members 
generally covered a well-defined area, of 5° to 10° 
diameter, whose center was chosen as radiant for 
the shower. Extraneous meteors within the radi- 
ant area were discarded when their velocity and 
time of occurrence indicated that they could not 
be shower members. For diffuse and relatively 
rapidly moving radiants, such as those of the 
Taurid streams, the data were not considered 


26 


sufficiently accurate to warrant the inclusion of 
straggling members. We have certainly under- 
estimated the strength of such streams. 

Table I compares our mean values with those 
accurately determined by Jacchia (Whipple 1954 
and unpublished) for the orbital elements, the 
velocity, and the radiant. For these last two we 
give the r.m.s. deviation from our mean value. 
This scatter should be interpreted as resulting 
from both the errors of measurement and the 
real spread in the shower parameters. The com- 
parisons show that our data adequately define 
the known showers and that the average orbital 
elements have at least statistical significance. 

After the known showers were identified, eight 
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groups of meteors remained that could be dis- 
tinguished as possible new streams. One of the 
eight, a high velocity stream occurring in mid- 
December, had already been discovered inde- 
pendently by Jacchia and will be discussed else- 
where. We deliberately made this initial search 
without knowing the location of theoretical or _ 
visual radiants to avoid prejudicing our choice. 
These new streams, although not spectacular, 
were in some cases more obvious than several of 
the known showers. We could accept the group 
of three meteors’in-Leo Minor only because this 
small, tight cluster appears apart from the ma- 
jority of the radiants. Such a group near the 
ecliptic, for example, would have passed un- 


TABLE I, METEOR SHOWER PARAMETERS DETERMINED BY TWO METHODS * 


Shower 1/ (17a) é q ” g 
Quadrantids 2.66 0.633 0.976 I71.2 282.4 
3.42 0.715 0.974 167.9. 282.2 
Perseids —29.4 1.03 0.950 153.6 138.9 
20.8 0.955 0.936 I51.2 138.1 
Orionids —59.9 1.01 0.580 80.6 28.8 
7.70 0.930 0.539 86.8 29.8 
S. Taurids 1.199 0.790 0.420 108.0 44.6 
2.30 0.835 0.380 III.9 45.1 
Geminids 1.41 0.900 0.142 323.9 260.8 
1.39 0.899 0.140 324.3 261.2 
a Capricornids 2.67 0.778 0.594 267.0 126.6 
2.36 0.756 0.577 270.2 124.1 
N. Taurids 2.42 0.851 0.364 292.2 227.7 
2.14 0.849 0.323 208.4 221.8 


* The first entry for each shower was determined by the approximat 
by the accurate reduction technique of Whipple and Jacchia is 


r.m.8. 
Radiant No.of r.m.s. Pee 
© Veo @ 8 cases Arad. lo 
93-6 (72.3 42.7 230 849 12 3°2 2.4 
90.I 73.8 42.7 230 «= 48 I 
292.5 112.7 60.8 44 #59 40 4.3 4-5 
289.3 113.7 60.4 46 = 58 II 
109.4 163.7 68.2 93 16 39 3.4 3.6 
116.5 163.2 66.5 94 16 2 
152.6 5-8 28.1 53 14 6 3.3 as 7. 
156.9 5-4 30.2 51 14 5 
224.7 23.2 36.4 TI25) 5 32 51 1.8 L.> 
225.6 24.0 36.5 113 32 19 
33-6 7.4 25.2 305 —I0 8 4-4 5.8 
34.3 7.1 25.1 304 —I0 6 
159.9 2.9 30.7 56 8622 10 2.6 4.6 
160.2 3-2 31.3 52 21 3 


method of McCrosky (1957); the second entry 


@ 


1957). 


TABLE II. PARAMETERS FOR 7 NEW PHOTOGRAPHIC SHOWERS AND 3 POSSIBLE COMETARY ASSOCIATIONS * 


62 cr No. r.m.8. 
te Radiant of rms. AV/V 
Shower 1/(iJa)- e q « 2 x s ofmax. duration Vio @ 8, cases Arad, of 
Coma 
Berenicieds —15.6 1.04 0.548 263.3 206.7 199.8 136.8 Jan. 17 Jan. 13-23 65.7 187 18 6 2 
— 5.99 1.08 0.4 270.5 209.8 210.3 130.6 64.7 185.5 20.5 z a4 ss 
Comet 19131 I 0.407 279.3 302.9 222.2 120.5 
- K Serpentids 8.04 0.944 0.449. 275.0 14.0 289.0 63.4 Apr. xr Apr. 1-6 46.4 231 18 4 3.3 4.2 
42.7 0.990 0.419 279.7 15.8 295.5 63.2 45.7 231.6 15.9 
@ Virginids 2.66 0.729 0.719 71.8 225.1 206,2 2.6 May 5 May 1-9 22.2 21s —12 9 5.8 9.0 
2.15 0.671 0.707 75.3 225.4 300.7 0.5 21.4 214.7 —14.8 I 
Leo Minorids 41.8 0.984 0.648 107.2 210.0 317.2 127.8 Oct. 2. Oct. 22-2 63.8 160 2 2.0 i) 
58.6 0.989 0.635 105.8 210.9 316.7 124.0 - 5 63.0 162.2 3s I : - 
Comet 1739 I 0.674 104.8 207.4 312.2 124.3 
e Geminids 14.1 0.944 90.790 234.9 205.6 80.4 171.4 Oct. 17 Oct. 17-27 70.1 102 26 6 3.5 1.6 
20.8 0.971 0.773 230.8 208.8 85.6 173.0 70.6 104.2 26.5 I 
uw Pegasids 4-35 0.778 0.965 199.4 229.7 60.1 8.3 Nov.1r Nov.11 16.1 340 23 
. s : fs : 5 5 a 3 2.3 6.7 
3.86 0.748 0.974 195.7 229.7 65. 8.1 15. 8 * 
Comet 1819IV 2.96 0.699 0.892 350.1 77.4 67. 9.1 oe er gre 
é Arietids 2.34 0.638 0.848 229.9 259.2 129.1 1.3 Dec. 8 Dec. 8-13 17.1 51 ar bi 3.8 4-7 


* When given, the second set of values for each shower is the result of accurate reductions by Jacchia (Whipple 1954 and unpublished). 
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noticed. Table II gives the mean orbital elements, 
velocities, and radiants for each of the new show- 
ers. Again, where possible, these data are com- 
pared with the elements of a member of the same 
shower which has been accurately reduced. 

Some uncertainty exists in the node for the 
a-Virginid stream. The members fall into two 
groups whose nodes differ by 180°; this is not 
surprising, however, in view .of the very small 
inclination of the orbit. 

We looked for cometary associations for the 
seven new radiants and definitely established 
three such links, a fraction comparable to the 
number of associations that have been made with 
the previously known photographic showers. 
The comparisons between the elements of shower 
and of comet appear in Table II. 

Our identification of the Coma Berenicieds 
with Comet 19131 receives support from the ten- 
tative connection suggested by Whipple (1954) 
between this comet and Harvard Meteor Number 
1918, the member of our shower whose elements 
were accurately reduced, as shown in Table II. 
The comet’s orbit is somewhat uncertain. David- 
son (1920) predicted a meteor stream associated 
with this comet at a = 187°, 6 = 22°, on Janu- 
ary 25. 

The association with Comet 1739 strengthens 
our selection of the Leo Minorids as a meteor 
stream. Denning (1925) had already suggested 
an association between this comet and a radiant 
at a = 155°, 6 = 40°, but had designated this 
association as ‘“‘probable or possible,’’ rather than 
“accurate or accepted.” The close correspond- 
ence between separate orbital elements in our 
data indicates an almost certain relation. 

Comet 1819IV is included in Porter’s (1952) 
list of comets approaching to within 0.1 a.u. of 
the earth. The predicted radiant for an associ- 
ated meteor stream is a = 330°, 6 = —39°, on 
January 9. A near identity of the longitude of 
perihelion, 7, exists between the » Pegasids and 
Comet 1819IV, but not between and w indi- 
vidually. Tisserand’s criterion for the identity of 
the orbits is satisfied. This short-period comet 
was subject to repeated perturbations by Jupi- 
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ter, so it is not unreasonable that the nodes of the 
orbit of the meteor stream should be the reverse 
of those of the original comet orbit, while Tis- 
serand’s criterion should be more applicable than 
for most meteor-comet identifications. 

The radiants listed by Denning (1899), Opik 
(1934), Hoffmeister (1948), McIntosh (1935), 
and Olivier (1921) were compared with the new 
radiants to reveal coincidences. The e-Geminid 
shower occurs, with varying degrees of certainty, 
in all the lists except those of McIntosh, whose 
observations were of southern radiants, and of 
Hoffmeister. If two closely associated groups of 
radiants in Denning’s list are combined, a reason- 
able argument can be made in favor of the ex- 
istence of the a Virginids. This shower is also 
present in the McIntosh data. The remaining 
five showers are either not found in the visual 
records or are so weak that their identity is ques- 
tionable. Their absence may be attributable, in 
part, to the large errors of visual observations. 
However, showers that show large yearly varia- 
tions in the number of members might be de- 
tected, by chance, in one year’s observations. 
If the observations covered a period of years, on 
the other hand, as do many of the visual cata- 
logues, the shower might be obscured by the 
accumulated sporadic background. 

The fact that nearly half of the new streams 
are associated with comets adds to our confidence 
in the results presented here. 


REFERENCES 


Davidson, M. 1920, M. N. 80, 739. 

Denning. R. M. 1899, Roy. Ast. Soc. Mem. 53, 203. 

——. 1925, Hutchinson's Splendour of the Heavens (New 
York: Robt. M. McBride and Co.), p. 431. 

Hoffmeister, C. 1948, Meteorstréme (Weimar: Verlog- 
Werdenund Wirken). 

McCrosky, R. E. 1957, Smithsonian Contr. to Astroph. 1, 


215. 

McIntosh, R. A. 1935, M. N. 95, 709. 

Olivier, C. P. 1921, Pub. McCormick Obs. 2, Part IV, 457. 

Opik, E. J. 1934, Cire. Harv. Coll. Obs. No. 388. 

Porter, J. G. 1952, Comets and Meteor Streams (New York: 
john iley and Sons), p. 92. 

Whipple, F. L. 1954, A. J. 59, 201. 

Whipple, F. L. and Jacchia, L. G. 1957, Smithsonian Contr. 
to Astroph. 1, 183. 


28 THE ASTRONOMICAL JOURNAL 


64, No. 1266 


COLOR-MAGNITUDE DIAGRAM FOR THE GLOBULAR CLUSTER M22 


By HALTON C. ARP anp WILLIAM G. MELBOURNE 


Mount Wilson and Palomar Observatories, Carnegie Institution of Washington, California Institute 
of Technology, Pasadena, California 


Received October 13, 1958 


Abstract. Photovisual and photographic magnitudes of more than 400 stars in the globular cluster M22 have been 
measured. The resulting color-magnitude diagram extends to slightly fainter than the RR Lyrae domain and is generally 
similar to those for other globular clusters, in particular most resembling the diagram for M13. 

Using the customary assumption that the mean magnitude of the RR Lyrae stars is Mp, = 0.0, an apparent distance 
modulus of 14.0 mag. and a reddening of 0.5 mag. are derived. The distance of the cluster is, therefore, 3.3 kpc and 


its brightest stars have M, = —3.0 and C.J. = +1.6 mag. 


1. Introduction. The globular cluster M22 is 
one of the largest and brightest in the sky and, 
appropriately, was one of the first whose stars 
were analyzed from the standpoint of their colors 
and magnitudes (Shapley 1930). That early color- 
magnitude array, however, was different from 
diagrams for any globular clusters known at the 
time or measured subsequently. As a possible 
explanation for this it was suggested that the 
M2z2 color-magnitude diagram was continuous 
with that of the rich Sagittarius star field which 
lies in the same direction as the cluster. 

The present investigation was undertaken, 
then, for the following reasons: 1) To measure 
the color-magnitude diagram of this nearby clus- 
ter on a photometric system which would enable 
it to be compared accurately to other clusters. 
2) To investigate the star field in the vicinity of 
M22, particularly with the view of obtaining 
information on the kind of stars occurring in this 
region of space which lies close to the galactic 
nucleus. 

2. Observational procedure. Photoelectrically 
measured standard stars are listed in Table I and 
identified in Plate I. The photoelectric measures 
were made in the same observing season as those 
in the Andromeda Nebula (Arp 1956) and the 
apparatus and techniques, therefore, were the 
same as described in that paper. 


TABLE I. PHOTOELECTRIC STANDARDS IN M22 


No. of 

Star ve B-V measures 
al 8.62 1.87 6 
a 10.55 39 6 
b 13.76 1.09 3 
c 14.42 1.63 2 
d 15.55 1.29 I 
e 15.14 30 2 
f 12.82 1.28 I 
g 12.22 1.41 I 
16.17 1.44 I 


Photographically measured stars are catalogued 
in Table II and identified in Plate II. Since the 
photographic‘plate and filter combinations used, 
103a-D + GGi1 and 103a-O + GGI, were on 
the m,., Mz, system, the photoelectric standards 
were converted by: 


“ty = B—0.27-+0.19 (B—V) 


MNy= V C.I. = My, —Mge- 


The remainder of the procedure is the same as 
described in an earlier paper on globular clusters 
(Arp 1955). 

3. The color-magnitude diagram. Figure 1 shows 
the color-magnitude diagram for M22. If the gap 
in the ‘‘horizontal sequence’ were well defined 
by stars both to the red and blue, as in M3 and 
Ms, there would be no difficulty in assigning a 
modulus and reddening to the cluster. In the 
absence of a bounded red side to the gap, it would 
still be well defined if the beginning of its blue 
end were approximately horizontal, as in M2, 
Mio and M15. 
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Figure 1. Color-magnitude diagram for M22. 
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Plate I. Photoelectric standards identified by letters. 
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TABLE II. PHOTOGRAPHIC MEASURES OF STARS IN M22 
(See Plate I and Plate II for identification) 


Star Vv CI Star 4 C.1. Star V C1. ‘Star Vv G.I Star Vv Cle 
e 15.13 .14 61 13.66 -93 Pal BECK ie SEO Sector III 70 14.41 eau 
f 123021 «25 62 12.29 -69 28 14.44 °&21.01 - Th. 28 5 71 14.95 89 
g 12:23 -1.40 63 14.02 OI 29 pT h27; .90 3 oF ae [Zan 04.28 22 
a 10.55 -22 64 14.48 .80 30 812.31 1.30 nor 0 2.19 73. 14.50 1.62 
b 13.76 1.03 O54. 798 Lull 31 12.05 1.43 3 : vf 1.05 homed Aste. 13 
c 14.42 1.67 66 14.17 1.03 3214.84 II 2 Te 81 75a eset 1.05 
d E5e5 5 eed 27, es athe .89 33 ve ne rs 12.85 1.17 3 ae I a 

12.69 1.11 34 14.82 aI : 

Sector I 69 14.43 84 35 14.80 .12 : as = Semel Sak Bes elie 
PPLigieT0= 260 70 14.92 .80 36 =-14.25 83 9 13.85 ‘99 49> 213.40) = ln? 
2 14.49 ak 71 15.13 1.07 37 14.94 -96 1013.58 192 80 §=14..38 73 
3. «18.11 “82 fae 214523) 31.50 38 14.88 1.44 II 13.32 2.01 81 13.18 .80 

(4 14.41 28 73 «=+13.72 -98 39 ~=15.01 B13 12. 11.55 1.77 82 13.72 1.01 

5 14.83 175 74 13.90 95 40 15.03 18 33 14.21 ‘41 33 14.45 15 
6 14.94 “15 75 12.02 1.24 41 14.48 .88 14 10.96 2.11 4 13°37; QI 
7 76 14.42 528 42 12.95 1.06 15 11.17 1.93 85. 13.71 94 
8 12.05 1.48 VPA 28) 21630 43 14.76 87 16 14.85 136 ee 12.61 as 
9 14.01 86 78 13.92 1.55 44 13.14 1.06 17 12.93 1.17 Wie 130223 ele 
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TABLE I1 (continued) 
Star V C., Star V Cts Star V Cal. Star V CI. Star Vv C.1. 
Sector IV A2e 14.48 -99 57. 14.10 -93 72) AL.35) 232 87- 13-32) 920.05 
28 14.79 oe AZ 2134530 9 1800 581 Ese Save he Ofna 7 smneel 2200 merst ale 3 er ood 
29 14.29 “94 44 15.05 -94 59) 2 Tl." 78 263 74 14.49 .86 2 ae a 
30. 14.47.03 45—«C«4: 79.4: OO G27 T0275 14.18 | 35 9 
31 «12.72 «1.28 46 15.01 ror 61. 12:71 4.32. 76 12.38 1.38 a es soe 
32 14.65 81 47 14.40 .89 62 13.34 1.14 77 14.74 -20 93 14.21 aaa 
33 12.62 1.29 48 14.33 -92 63 14.30 -99 78 14.85 1.05 904 14.60 .30 
34) 34°38.) 61..08 49 13-51 -99 64) T1408") 1.02 FOS T4ALSO! Walz 95 > 13.66 1.50 
35. 14.92 -98 50}; 314520 -38 65 13.86 . 83 80 14.37 -66 96 13.24 1.09 
36> ..14.86> 1500 51 14.21 ay) 66 13.24 1.06 81 14.90 .82 97 II.20 2.04 
37. 14.26 .93 52. 5) 14220 .97 G7 2225557 1.23 82 14.00 .88 98 14.60 27 
38 «14.88 -95 53 14.01 -40 63). 5357. OI 83 14-31 .OI 99 ‘Igte@guee srs 
BOmnr ise 20 ok sOz 54 14.43 -99 69 14.65 .19 84 14.20 . 1256 : 100%) 1268s est 
40 13.08 1.15 55 14.93 -QI 70 13.94 -93 85 13.55 1.02 Io! 13.03 1.08 
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The blue sequence in M22, however, shows 
considerable slope. This is true of only one other 
globular cluster measured previously, M13. The 
blue sequence in Figure 1 has therefore been 
matched to the corresponding feature in M13, 
and the gap in Figure 1 has been drawn to corre- 
spond to the position of the RR Lyrae stars 
measured in M13. 

This procedure yields a modulus of m — M = 
14.0 and a reddening of Ec.r. = 0.5 mag. If we 
take the other extreme case, i.e., that the gap be 
placed exactly adjacent to the end of the blue 
sequence, we derive m — M = 14.2 and Ec.z, = 
0.4.mag. Although the first set of values has been 
adopted here as most probable, the possible 
values allowable could range to the latter limit. 

In a diagram such as M22, the best method of 
exactly placing the zero point is, of course, to 
measure the RR Lyrae Cepheids themselves. 
Since the cluster does contain a number of RR 
Lyrae Cepheids the modulus quoted here, though 
presumably accurate to within 0.1 -mag., is ca- 
pable of being improved. “ 

The general form of the color-magnitude dia- 
gram for M22 is obviously that of a normal 
globular cluster. Placing the zero-point as above, 
yields M, = —3.0and C.J. = +1.6 mag. for the 
brightest stars in the giant branch. The sequences, 
‘then, lie within the range found for other globu- 
lar clusters and the first array, discussed in the 
introduction, must be regarded as due to a scale 
error in the magnitudes. 

It is encouraging to note that when the color- 
magnitude diagram for M22 is matched in the 


region of the blue sequence with that of M13, 
the giant branch of M22 then coincides exactly 
with the giant branch of M13. The two clusters 
seem to be, on the whole, closely the same except 
_for the following small differences: 1) The giant 
branch of M22 comes down to M, = 0.0 at a 
slightly bluer color index than in M13 and then 
breaks more sharply into the vertical sub-giant 
branch. The shape of this feature most resembles 
the giant-sub-giant transition phase in M3, al- 
though the whole feature is moved more thar 
0.2 mag. bluer in color index than in M3. 2) The 
well-defined bifurcation in the giant sequence of 
M13 is not present in the M22 giant seque.ice. 
In fact, the bifurcation in the M22 giant branch 
is less noticeable than in any of the other seven 
globular clusters previously measured on a simi- 
lar basis. 

In agreement with the general similarity be- 
tween M22 and M13, M22 then represents an- 
other case in which the parallel displacement in 
color index of the giant branches, observed in M3 
and Mo2, and in others, is violated. 

The spectra of the brighter stars in M22 are 
described by A. J. Deutsch (unpublished) as not 
so excessively metal-poor as those in M15 and 
Moz2, but similar instead to stars in clusters such 
as M13 and M3. 
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Abstract. Analysis of the colors and magnitudes of 300 stars in a direction 15’ north of M22 is made. The results indi- 
cate that within this volume of space, which passes within 2.3 kpc of the galactic center, the majority of redder stars 
are less luminous than globular-cluster giants of the same color index. These stars appear to be nearer in luminosity to 


old galactic-cluster giants, 


With the 60-inch telescope diaphragmed to 32 
inches, photovisual and photographic plates 
(103a-D + GGi1 and 103a-O + GG13) were 
centered on an area 15’ north of the globular 
cluster M22. Since the photoelectric standards 
listed in the preceding paper were within the 
photometric field of these plates, colors and mag- 
nitudes of more than 300 stars in an area 10’ by 
15’ could be measured relative to those standards. 
The resulting color-magnitude diagram is shown 
in Figure 1. 

It can be readily seen, by noting the position 
of the schematic diagram for M22, that the field 
stars in the vicinity have no appreciable effect 
on the color-magnitude diagram for M22. The 
results obtained for M22 in the preceding paper, 
however, enable a rudimentary analysis of the 
kind of stars which occur in this sample volume 
of space directed towards the galactic center. 

Since the plate and filter. combinations used in 
this investigation give magnitudes on the B, V 


Figure 1. Color-magnitude diagram of field 15’ north of 
M2z2. Upper solid line is schematic color-magnitude dia- 
gram of M22; lower dashed line is how M22 would appear 
at distance of nucleus with the same reddening. 


system, the only data which need be changed 
from the preceding paper are those dealing with 
the color indices. We have for the globular clus- 
ter, M22: 


m— M = 14.0 1 = 337° b= —9° 
Es_v = Ec.r. X 0.81 = 0.4 mag. 


Using the relation between reddening and ab- 
sorption (Blanco 1956) we can compute the ab- 
sorption, the distance of the cluster from the sun, 
and its height above the plane: 


A, =1.4mag. 7 = 3.3kpc 2 = 500 psc. 


Taking the distance to the galactic center as 8.7 
kpc (Baade 1950), it is calculated that the vol- 
ume sampled in the field under investigation 
passes 2.3 kpc from the galactic center. The 
geometry of the position of M22 relative to the 
center is shown in Figure 2. 

The place in Figure 1 which would be occupied 
by type II (globular cluster) stars at the distance 
of the nucleus may be computed by assuming the 
following about the reddening: 1) The absorption 
in the direction of M22 lies entirely between 
M22 and the sun, i.e., the obscuring material 
does not extend more than 500 parsecs above the 
plane of the galaxy. 2) The absorption across the 
face of the star field investigated is the same as 
across M22. This gives a total obscuration of 
A, = 1.4 mag. to the center. Two more assump- 
tions are made: 3) The distance to the nucleus 


NUCLEUS 


Figure 2. Position of M22 and background star field. 
Distances are in kiloparsecs. 
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is 8.7 kpc and 4) the absolute magnitude of the 
RR Lyrae stars is M,, = 0.0 mag. 

It is then possible to compute how the sche- 
matic diagram of a globular cluster at the dis- 
tance of the nucleus looks, and this is shown in 
Figure 1 by the dashed line. It is apparent. that 
the stars measured in this region of space do not 
show any preference for the type II color-magni- 
tude diagram. 

The stars which appear in Figure 1 really do 
reside in the galactic center regions as can be 
seen from two arguments. 1) Considering only 
the stars redder than B— V = 1.5 mag., it can 
be computed that nearby main-sequence stars, 
even if unobscured, would be bright enough to 
appear in Figure 1 only if they were closer than 
about 200 parsecs. On the basis of star density 
near the sun (Bok 1937 and Allen 1955) no more 
than about five such stars would be predicted to 


occur in the present field. The 95 which do ap- - 


pear brighter than V = 16 mag. and B—V = 
1.5 in Figure 1, must then be giants at large 
distances from the sun. 2) Independently, Seares’ 
(1925) average counts at this latitude would pre- 
dict about 20 stars of all color indices brighter 
than V = 16 mag. Finding so many more in the 
longitude of this field again indicates the stars in 
Figure I originate well beyond the solar neigh- 
borhood and are giants contributed by the den- 
sity build-up toward the galactic center. 

Figure 1 then demonstrates that the stars in- 
vestigated here, which lie roughly within several 
kiloparsecs of the center, are less luminous than 
globular-cluster giants. The luminosity of these 
stars is more appropriate for old galactic-cluster 
giants. 

If the absolute magnitude of the RR Lyrae 


stars is fainter than M,, = 0.0 mag., which is. 


probable (Arp 1958), or if the nucleus is some- 
what more distant, or if the reddening to it in 
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this field is somewhat more than computed here 
then the dashed globular-cluster schematic wil 
be moved downward slightly. Then some of the 
stars observed here could be reasonably identi. 
fied as globular-cluster giants at the distance o! 
the nucleus. The majority of the stars woulc 
still be under-luminous, however, and closer tc 
galactic-cluster giants. This result is hard to rec: 
oncile with Baade’s (1950) high density of RR 
Lyrae variables in the nucleus. It should be 
noted, however, that Baade’s observed densitie: 
indicate that the type II population is very con: 
centrated to the center. To illustrate, 80 per cen’ 
of his RR Lyrae variables lie within 2 kpc of the 
center. The volume of space sampled here passe: 
further than 2 kpc from the nucleus. It may be 
that the old type I, the old galactic cluster popu 
lation, is actually dominant in the region investi 
gated here. 

The present-day concept, that the type II 
globular cluster population, has a ‘“‘halo’’-like 
distribution, however, is in contradiction wit 
Baade’s extreme concentration of the RR Lyrae’ 
to the center and makes the above interpretatior 
uncertain. It is apparent that further observa 
tions are needed to understand the compositior 
of the central regions of our galaxy. Color-mag 
nitude diagrams, particularly in Baade’s NGC 
6522 region, and three-color work in these anc 
other nuclear regions are required. 
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The Astronomers’ Universe, by Bart J. Bok, published by New York: Cambridge University Press, 
32 East 57th St., New York 22, N. Y. Price $3.75. 

The Threshold of Space, edited by M. Zelikoff, published by Pergamon Press Inc., 122 East 55th St., 
New York 22, N. Y. Price $15.00. 

L’Exploration des Galaxies Voisines, par les méthodes optiques et radio-électriques, by Gérard de 
Vaucouleurs, published by Masson et Ci* Ed., Paris, France. Price 1.600 fr. 


NOTICES 


The American Institute of Physics has undertaken the complete translation of all issues of the 
“Astronomical Journal of the Academy of Sciences of the USSR,”’ commencing with the 1957 issues. 
As the title suggests, this journal covers various problems in astronomy, including solar activity, 
stellar studies, spectroscopic investigations and radio astronomy. 

There are six issues annually; subscription price: $25.00 domestic, $27.00 foreign ($10.00 and 
$12.00 respectively for libraries of nonprofit degree-granting institutions). Subscriptions and in- 
quiries should be addressed to: Translation Journals—American Institute of Physics, 335 East 45 
Street, New York 17, N. Y. 


A Symposium on Boundary Problems in Differential Equations, with special reference to recent 
developments in this field, will be held by the Mathematics Research Center, United States Army, 
at the University of Wisconsin April 20-22, 1959. Invited speakers, about twenty in number, will 
each present a thirty-minute paper. Both ordinary and partial differential equations will be con- 
sidered, the emphasis to be upon methods that are potentially adapted to computation. Between 
lectures there will be discussion periods and coffee breaks. The proceedings of the Symposium will 
be published. 

Among the speakers, European and American, will be 


L. Collatz R. Bellman K. Friedrichs 

G. Fichera G. Birkhoff P. Garabedian 
L. Fox H. Bueckner B. A. Troesch 
W. T. Koiter R. Courant R. Varga 

J. Schréder J. B. Diaz C. Wilcox 

I. N. Sneddon J. Douglas D. Young 


Persons interested in attending the symposium may receive the program and other details by 
writing to 
R..E. Langer, Director 
Mathematics Research Center, U. S. Army 
1118 West Johnson Street 
Madison 6, Wisconsin 
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EDITORIALS 


REDUCED SUBSCRIPTION RATE FOR MEMBERS OF 
NATIONAL ASTRONOMICAL SOCIETIES 


By action of the Council of the American Astronomical Society, bona fide resident members of 
the various national astronomical societies (e.g. Fellows of the Royal Astronomical Society resident 
in the United Kingdom and Fellows of the Royal Astronomical Society of Canada resident in Canada) 
will be allowed the opportunity of subscribing to the Astronomical Journal at the rates currently 
available to members of the American Astronomical Society. Those who wish to take advantage of — 
these reduced rates, should write to the Editors of the Astronomical Journal stating that they are 
duly elected members of their respective national astronomical societies and that they wish to sub- 
scribe to the Astronomical Journal at the reduced membership rates. 


PUBLICATION OF PAPERS IN RADIO ASTRONOMY 


The accelerated development of radio astronomy in America and the need for rapid publication — 
in this field has raised doubts in the minds of some radio astronomers whether the existing astro-- 
nomical journals in this country are adequate for their needs. For this reason, we should like to © 
state that the Astronomical Journal heartily welcomes authors to submit papers on radio astronomy — 
and assures them that their papers will be given every consideration. This has always been the 
policy of the Astronomical Journal: during 1957 and 1958, the Astronomical Journal published several f 
papers on radio astronomy ; and several are in course of publication. 

Radio astronomy is a rapidly developing branch of astronomy ; and any schism between thia’ 
newer branch and the older branch of optical astronomy is to be greatly regretted. For this reason, © 
the Astronomical Journal, already established as one of the principal astronomical journals, would ' 
like to see that papers on radio astronomy, commensurate with its growing importance, will be © 
represented in the Astronomical Journal. 

Dirk BROUWER 
HARLAN J. SMITH 
* Editors 
Yale University Observatory 


